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ABSTRACT 

We present initial results from the low-latitude Galactic plane region of the High Time Res¬ 
olution Universe pulsar survey conducted at the Parkes 64-m radio telescope. We discuss the 
computational challenges arising from the processing of the terabyte-sized survey data. Two 
new radio interference mitigation techniques are introduced, as well as a partially-coherent 
segmented acceleration search algorithm which aims to increase our chances of discovering 
highly-relativistic short-orbit binary systems, covering a parameter space including potential 
pulsar-black hole binaries. We show that under a constant acceleration approximation, a ratio 
of data length over orbital period of sa 0.1 results in the highest effectiveness for this search 
algorithm. From the 50 per cent of data processed thus far, we have re-detected 435 previously 
known pulsars and discovered a further 60 pulsars, two of which are fast-spinning pulsars with 
periods less than 30 ms. PSR J1101—6424 is a millisecond pulsar whose heavy white dwarf 
(WD) companion and short spin period of 5.1 ms indicate a rare example of full-recycling via 
Case A Roche lobe overflow. PSR J1757—27 appears to be an isolated recycled pulsar with a 
relatively long spin period of 17 ms. In addition, PSR J1244—6359 is a mildly-recycled binary 
system with a heavy WD companion, PSR J1755—25 has a significant orbital eccentricity of 
0.09, and PSR J1759—24 is likely to be a long-orbit eclipsing binary with orbital period of 
the order of tens of years. Comparison of our newly-discovered pulsar sample to the known 
population suggests that they belong to an older population. Furthermore, we demonstrate that 
our current pulsar detection yield is as expected from population synthesis. 

Key words: stars; neutron — pulsars; general — methods; data analysis — surveys. 


1 INTRODUCTION 

Pulsars are extraordinary natural laboratories with extremely high 
density and gravity impossible to re-create on Earth, hence they 
provide exclusive insights to a rich variety of fundamental physics 
and astronomy. To discover more pulsars we have performed the 
High Time Resolution Universe (HTRU) survey; a blind survey of 
the southern sky w ith the 64-m Parkes radio telescope in Australia 
dKeith et al.ll2010h and a twin survey of the northern sky with the 


100-m Effelsberg radio telescope in Germany (the HTRU-North 
survey; iBarr et aOllOlSh . The HTRU survey uses multi-beam re¬ 
ceivers and backends constructed with state-of-the-art technology, 
providing unprecedented high time and frequency resolution, al¬ 
lowing us to probe deeper into the Galaxy than previous efforts 
with these telescopes. We split the survey into three regions of the 
sky, namely the low-latitude Galactic plane survey, the medium- 
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latitude survey, and the high-latitude all sky survey, each tailored to 
achieve specihc science goals. 

The low-latitude Galactic plane regio n is where the most rel - 
ativistic binaries are expected to be found jBelczvnski et al.|[2002h . 
Pulsars in tight binaries orbiting other compact objects, for ex¬ 
ample neutron stars and, potentially, black holes, are of great in¬ 
terest as their strong gr avitational fi elds provide the best tests of 
General Relativity (GR: IWexll2014h and other theories of gravity 
dpreire et alj2012h . The best example of such a binary system so fa r 
is the double pulsar system teurgav et al.ll2003l : lLvne et alj|2004ll . 
The double pulsar has been used to obtain four independent tests 
of GR and GR has passed the most stringent test, regarding the 
shape of the S hapiro delay, with a m easurement uncertainty of only 
0.05 per cent jRramer et al.ll2006^ . The number and the precision 
of GR tests increase as the binary systems to be discovered be¬ 
come more relativistic. Hence, one of the main aims of the HTRU 
Galactic plane survey is precisely the discovery and study of ultra¬ 
compact relativistic binary systems in short orbits. 

In this region of sky within Galactic latitude ±3.5°, we em¬ 
ploy the longest HTRU integrations of 72 min per pointing to 
maximise our sensitivity. This means that the HTRU low-latitude 
Galactic plane data set will be capable of revealing many pulsars 
that were not luminous enough to be detected by previous surveys. 
These pulsars might increase the sample of sources that glitch, 
whic h could lead to improv ed knowledge of the interior of neutron 
stars tespinoza et al.ll201lh . They will give us useful insights into 
the interstellar medium via the study of their dispersion and rota¬ 
tion measures to reveal a picture of the magnetic field of the Milky 
Way (see e.g.. lNoutsos et al]|2008l) . They also provide us with an 
important handle on the lower end of the luminosity distribution 
function of the Galactic plane pulsar population, valuable knowl¬ 
edge for the planning of survey strate^s for the next generation 
of radio telescopes, such as MeerKAlU and the Square Kilome¬ 
tre Arra)0(SKA). Furthermore, the archive of the HTRU Galactic 
plane survey will continue to produce science through future data 
re-examination. The long dwell time of these observations is also 
favourable for the detection of transient and nulling sources deep 
withi n the Galactic plane (see e.g., related work by IPetroff et al.l 

[20Ti) . 

However, the sheer volume of the high-resolution HTRU 
Galactic plane data set poses great challenges in data manipula¬ 
tion and analysis. Normally, long integration length (fint) provides 
an increase in sensitivity. Nonetheless, this is not exactly the case 
when it comes to searching for tight-orbit relativistic pulsar bi¬ 
naries: a periodicity search in the Fourier domain is the standard 
method employed in most pulsar surveys. However, the high orbital 
acceleration attained by fast relativistic binaries results in a Doppler 
shift in the spin frequency of the pulsar as a function of the orbital 
phase. The pulsar signal is thus smeared across neighbouring spec¬ 
tral bins of the Fourier power spectrum, hence a reduction in the 
sensitivity of the periodicity search. Furthermore, the width of the 
Fourier spectral bin is defined by 1 /tint. As a result, the longer the 
integration time, the larger the portion of the orbit we cover in a 
particular observation, but also the narrower the Fourier spectral 
bin. These two effects combined lead to more severe consequences 
for Fourier spectral smearing. For these reasons, the search algo¬ 
rithm of the Galactic plane survey presented in this paper requires 
significant modifications over that used in the medium- and high- 


^ http ://www. ska.ac. za/meerkat/index .php 
^ https;//www.skatelescope.org/ 


Table 1. Characteristic minimum detectable flux density (S 1400 min) for 
the HTRU Galactic plane survey. Considering normal pulsars and MSPs as 
two separate groups, we note the minimum, mean, and maximum of the 
duty cycle <5 of these two groups respectively, and we derive the respective 

51400 ,min- 


<5(%) 

MSPs 

‘S'l400,min (^TlJy) 

Normal pulsars 

5 (%) 5'i4QQ^min (niJy) 

min 0.28 

0.013 

0.014 

0.0030 

mean 11.54 

0.092 

4.21 

0.053 

max 65.31 

0.35 

57.29 

0.29 


latitude survey jKeith et al.ll2010h . Various search techniques tar¬ 
geting binary pulsars exist, which compensate for, or in some cases 
fully rec over, the reduction in det ectability due to orbital motion 
(see e.g.. lLorimer & Krameill200^ . Here we present an innovative 
segmented search technique which aims to increase our chances of 
discoveries of highly-accelerated relativistic short-orbit binary sys¬ 
tems, including potential pulsar-black hole binaries. We stress that 
the depth of the parameter space to which the survey data can be 
explored is highly dependent on the available computing resources. 
Optimisation of pulsar searching algorithms is thus crucial in the 
era of data intensive astronomy, and the computational challenges 
faced by the HTRU Galactic plane survey will be applicable for the 
planning of the SKA. 

In this paper we focus on the data processing and initial results 
of the HTRU Galactic plane survey. In Section |2] we describe the 
survey strategy and the observation set-up. In Sectionl^we present 
the data processing algorithm with emphasis on the newly incorpo¬ 
rated elements, namely two radio-frequency interference (RFI) mit¬ 
igation schemes (Section[3T| and an innovative ‘partially-coherent 
segmented acceleration search’ technique (Section[±2ll. We report 
on the re-detections of known pulsars in Section|4]and the discov¬ 
eries of 60 pulsars in Section|3 We highlight the individual pulsars 
of interest in Section[^and we compare the newly-discovered pul¬ 
sars with the known pulsar population in Section|7] In Section[^we 
compare our discovery rate with the estimated survey yield and in 
Sectionl^we present our conclusions. 


2 GALACTIC-PLANE SURVEY STRATEGY 

The Galactic plane survey covers a strip along the inner Galac¬ 
tic plane, with the central beam of all scheduled pointings be¬ 
tween Galactic longitude —80° < / < 30° and latitude |6| 
< 3.5° (see Fig. [Hi. The observ ational set-up is the same as that 
described in [Keith et ^ (l201(lh . To summarise, the HTRU sur- 
vey observations were mad e using the 20-cm multibeam receiver 
dStavelev-Smith et^ll996ll on the 64-m Parkes radio telescope. 
This receiver is designed for efficient sky surveying, allowing si¬ 
multaneous observations with its 13 receiver beams each separated 
by ~30' (Fig. [ID. Each receiver beam has a full width at half¬ 
maximum (FWHM) of 14'.4, and thus requires four interleaving 
pointings to cover an area in mosaic style (see Fig. [D. A band¬ 
width of 400 MHz is recorded for each beam, but the introduction 
of low-pass hardware filters against RFI generated by the Thuraya 3 
geostationary communications satellite means that we have a final 
usable bandwidth of 340 MHz centred at 1352 MHz. 

We employ the Berkeley-Parkes-Swinburne Recorder (BPSR) 
digital backend and we sample 1024 channels across the observ¬ 
ing bandwidth giving an intra-channel frequency resolution of 
0.39 MHz, at a sampling rate of fsamp = 64 p.s. These are signifi- 
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cant improvements from the previously most suc cessful pulsar sur¬ 
vey, t he Parkes Multibeam Pulsar Survey fPMPS: lManchester et alj 
l200lh . where the frequency resolution was 3 MHz and time resolu¬ 
tion 250 ps. These mean that the HTRU survey has high frequency 
resolution to help the removal of interstellar dispersion and is sen¬ 
sitive to the transient sky down to 64 p.s. In addition, the observa¬ 
tions were decimated to 2 bits per sample. This is another advan¬ 
tage compared to PMPS where its 1-bit digitisation might have re¬ 
duced sensitivity for weak signals and bright radio bursts. Based on 
these survey specifications, we can use the radiometer equation (see 
Equation ([5]l) to calculate the characteristic minimum detectable 
flux density at 1.4 GHz ( 51400 , 111111 ) for the HTRU Galactic plane 
survey. This quantity depends on pulsar duty cycle (5), hence we 
have derived the characteristic 5 i4oo,min with respect to the mini¬ 
mum, mean and maximum 5 of al l published pulsar data from the 
ATNF Pulsar CataloguJ^(PSRCAT: lManchester et alj2005h . We list 
in Table[T]the corresponding 5 i4oo,min distinguishing between nor¬ 
mal pulsars and separately for millisecond pulsarfl (MSPs). 

Observations for the HTRU Galactic plane survey at Parkes 
took place between November 2008 and December 2013, which 
comprise 1230 scheduled pointings each with 72-min long obser¬ 
vations with just over 2^® samples. The observations have their 
two polarizations summed, before being written to tape for stor¬ 
age. Several corrupted pointings (due to severe RFI contamination 
or hardware issues) were re-observed, hence finally 1246 pointings 
were recorded, resulting in 263 terabytes of observational raw data. 

3 PROCESSING 

Two processing pipelines are applied to the data. Our initial, 
‘standard search’ pipel ine, also known a s the HITRUN processing 
pipeline, is outlined in iKeith et al.l (l2010h . This algorithm follows 
th e typical procedures of pu lsar searching such as that described 
in lUorimer & Rramej j2005h . Firstly, spurious signals in the data, 
such as those created by radio frequency interference (RFI), are 
identified and excised. Next, the observation is dedispersed to com¬ 
pensate for the frequency-dependent delay caused by the free elec¬ 
trons along the line of sight. As the amount of dispersion is de¬ 
pendent on the a priori unknown distance of the pulsar to be dis¬ 
covered, we trial a wide range of potential dispersion measures 
(DMs) between 0 and 3000 cm~® pc, which sums to a total of 1069 
DM trials per data set. Each of the dedispersed time series is then 
Fourier transformed. We sum the second, fourth, eighth and the 
16th harmonic Fourier spectr a respectively (see ‘incoh erent har¬ 
monic summing’ technique in iTavlor & Hugueninlll96lj) and iden¬ 
tify significant signals creat ed by any pulsating sourc es. Based on 
the false-alarm probability dUorimer & Kramenl2005l) . this survey 
has a signal-to-noise (S/N) threshold of ~ 9. We have nonetheless 
inspected every potential pulsar candidate with a S/N above 8 by 
eye. A second pipeline, the ‘partially-coherent segmented acceler¬ 
ation search’, aims to improve the detectability of binary pulsars 
and is detailed in Section lT^ In addition, advancements regarding 
RFI mitigation are implemented and ar e presented in Sectionl3.ll 
A single-pulse related analysis (see e.g.. lBurke-Spolaor et al.l201lh 
has not yet been carried out and will be conducted in future data re¬ 
processing. 

To date, 615 pointings of the Galactic plane survey have been 
® http://www.atnf .csiro.au/people/pulsar/psrcat/ 

^ iLee et al] 1201 2h derived an empirical definition to classify MSP. For sim¬ 
plicity, we adopt a definition of P Si 30 ms for MSP throughout this paper. 



Figure 2. Histogram showing the statistics of the percentage of time sam¬ 
ples removed per observation, as a result of the automatically generated 
RFI mask. The blue line shows the conesponding cumulative distribution 
function. 

(Statistics taken from 414 observations, i.e. 34 per cent of the survey) 

processed, which is 50 per cent of the survey. The spatial distribu¬ 
tion of the processed data can be seen in Fig. [T] These processed 
pointings are not contiguous and only depend on availability of data 
at the location of the computing facilities. 


3.1 RFI mitigation 

iKeith et alJ ll201(]h described two RFI removal procedures em¬ 
ployed as part of the HITRUN pipeline: the removal of RFI-affected 
spectral channels in the frequency domain targeting narrow chan¬ 
nel interference, and the replacement of time samples contaminated 
by impulsive RFI with noise generated from random sampling of 
the uncontaminated surrounding data. Two extensions have been 
incorporated in the Galactic plane survey, both exploiting the fact 
that RFI is terrestrial hence usually non-dispersed (i.e., most promi¬ 
nent at DM = 0 cm~^ pc) and often appears in multiple receiver 
beams, in contrast to celestial sources which are point-like and tend 
to show up in only one beam, unless they are very bright. 


3.1.1 Time doma in 

Since November 2010 an automated scheme to generate empiri¬ 
cal time-domain RFI masks has been incorporated into the BPSR 
backend. T his RFI mitigation algorithm has been developed from 
the work of IKocz et al.l (l2012h . During each HTRU survey obser¬ 
vation, non-dispersed (i.e., at DM = Ocm~®pc) time series are 
output every 10 s from BPSR. Each of these 10-s time segments 
are auto- and then cross-correlated to form a covariance matrix and 
are subsequently decomposed into eigenvalues. We apply a thresh¬ 
old to identify RFI-affected time samples, namely a cut at 6a for 
any signal that appears in more than four beams and a lower cut at 
Aa if the signal is present in all 13 beams. These potentially con¬ 
taminated time samples are recorded to a log file, which is stored 
together with the un-corrected observation. The RFI mask is ap¬ 
plied to all beams of a single pointing during subsequent off-line 
data processing to replace the bad time samples with random noise. 
Fig. [2] shows the statistics of the percentage of masked time sam¬ 
ples per observation for the HTRU Galactic plane survey, which 
also serves as a measure of the quality of the survey data. Typi¬ 
cally, for each observation, ~ 0.03 per cent of the time samples 
are flagged as RFI-affected, and no observation has > 2 per cent 
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Figure 1. The spatial distribution of the processed observations from the HTRU Galactic plane survey. Grey circles denote the 1230 scheduled pointings for 
the complete Galactic plane survey. Blue circles denote pointings which have been processed with the ‘standard search’ pipeline (amounting to ~14 per cent 
of all pointings), whereas red circles ai'e pointings which have been processed with the ‘partially-coherent segmented acceleration search’ pipeline (amounting 
to ~35 per cent of all pointings). A small portion of 156 beams of observations have been processed with both pipelines to check for compatibility between 
the two algorithms and are represented by purple circles. Yellow crosses show the locations of previously known pulsars and green crosses show the 435 
re-detected pulsars. The black stars show the 60 pulsars discovered thus fai'. 
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of masked time samples. For the observations taken before this im- presented here, which from our retrospective checks show that they 

plementation, we apply the tim e domain RFI miti gation technique would not have been found otherwise, 

as described in Section 4.1.1 of iKeith et al.l (l2010h . 


3.1.2 Fourier domain 

The presence of any periodic RFI is most readily identified in the 
Fourier domain. Nonetheless, great caution must be taken when 
flagging periodic RFI, to prevent accidentally masking a genuine 
pulsar signal which is also periodic. As part of the HITRUN pipeline, 
fixed periodic RFI lists have been applied to mask known RFI pe¬ 
riodicities, such as that of the Australian mains power supply at 
~50Hz and its harmonics. However, the RFI environment is tem¬ 
porally varying and a fixed list is far from ideal. We have investi¬ 
gated a more robust method of empirically identifying birdie pe¬ 
riodicities. For each observation, we create the power spectrum of 
each of the 13 beams from their time series at DM = 0 cm“® pc. 
A power threshold (Pthres) based on false-al arm probability can 
be cal culated with the following equation from iLorimer & Kramej 
(l2005h . 

Pthres = - In ( —i - ) , (1) 

y ^n^samp / 

where nsamp is the number of time samples used to create the 
power spectrum and for this survey is 2^®, which gives a Pthres ~ 
19. We then compare the power spectrum Fourier bin by Fourier 
bin, identifying Fourier frequencies which exceed "Pthres in more 
than four beams and flagging them as periodic RFI. Under this 
scheme, typically < 0.01 per cent of the Fourier spectrum is re¬ 
moved for each data set. This compares very favourably to the 
fixed periodic RFI list incorporated as part of the HITRUN pipeline, 
which contains Fourier frequencies associated with known RFI and 
always results in the removal of 0.57 per cent of the Fourier spec¬ 
trum. Note that for the acceleration search algorithm as described 
in Section as the time re-sampling at every acceleration trial 
would shift any RFI periodicities in the time series, one periodic 
RFI list per acceleration trial has to be created. 

Fig. 0 compares the effectiveness in mitigating periodic RFI 
by applying varying stages of the above-mentioned techniques on 
one of the survey observations. In all cases, the original uncleaned 
spectrum is plotted in red, and the spikes corresponding to harmon¬ 
ics of the Australian mains power supply at ~ 150 Hz can be clearly 
seen. Plotted in black are the spectra after each RFI cleaning pro¬ 
cedure, and in theory no spike should remain, as there is no pulsar 
in this observation. In panel (a) only the time domain mask as de¬ 
scribed in Section [3. 1.1 1 has been applied, which might have been 
effective for removing impulsive RFI but remains insensitive in the 
Fourier space. In panel (b) also the frequency channels with exces¬ 
sive po wer are removed as described in Section 4.1.1 of ! Keith et al.l 
( l2010h . but the result is almost identical to the previous panel. In 
panel (c) we apply the fixed periodic RFI list incorporated as part 
of the HITRUN pipeline, which turns out to be masking a part of the 
spectrum that is relatively RFI-free, but is not able to identify some 
of the neighbouring narrow spikes. Finally, panel (d) shows the 
result of applying the empirically-determined multi-beam Fourier 
domain RFI mitigation method presented here. Occasionally a few 
RFI-related spikes (for example those at ~ 148.8 and 148.9 Hz), al¬ 
though appearing significant by eye inspection, are still missed out 
by this RFI mitigation technique, as they are present in less than our 
chosen conservative threshold of four beams. Nonetheless, most of 
the prominent spikes have been successfully identified and masked. 
Indeed, this technique has enabled several of the pulsar discoveries 


3.2 Acceleration search 


The challenges related to the searches for binary pulsars have been 
mentioned in Section[T] To fully account for the effects of unknown 
orbits on pulsars in our data is an almost impossible task, as theoret¬ 
ically one should search five of the Keplerian orbital parameter^ 
as well as all possible DMs and spin periods, resulting in a large 
7-dimensional parameter space which would be extremely compu¬ 
tationally expensive. 

A more manageable strategy is to approximate any unknown 
orbital motion as a simple line-of-sight constant acceleration, i.e., 
v{t) = at, combined with a ‘time d omain resampling’ technique 
(see e.g.. |johnston & Kulkarniiri99lh . By quadratically stretching 
or compressing a time series by the amount dictated by a particular 
acceleration, the time series is re-binned into equal time steps in the 
rest frame of an inertial observer with respect to the pulsar in binary 
orbit. This resampled time series can then be Fourier transformed to 
coherently search for peaks in the power spectrum just like the stan¬ 
dard periodicity search. This ‘time domain resampling’ technique 
has been a frequent choice for previous pulsar surveys targeting 
bi nary pulsars. Notabl e examples are the 47 Tucana e observations 
bvIcWilo et al.l ( l2000h and the PMPS re-analysis bv lEatough et al.l 
( 12013 ). In the following we detail the implementation of such an 
acceleration search for the HTRU Galactic plane survey. 


3.2.1 The ratio of data length over orbital period, Torb 

The constant acceleration approximation is equivalent to the best¬ 
fitting tangent to a quadratic v{t) curve, and its effectiveness thus 
depends on the ratio of the integration length (fint) to the orbital 
period of the pulsar (Porb). Here we define this ratio to be Torb, 
where 

lint fsamp N tlFFT 

‘Kit- ™ 

Note that tint is the product of the time sampling rate, fsamp, and 
the number of samples used in the Fast Fourier Transform (FFT), 
jifft. This number should be a power of two for maximum com¬ 
putational efficiency of the FFTs. 

To determine the effectiveness of this constant acceleration am 
proximation versus varying Torb, we employed two observation|j 
of the double pulsar system PSR 10737—3039A as test data sets 
(refer to Appendix for details of this study). Fig. |4]is a plot 
of the S/N of four selected orbital phases from this study across 
varying Torb. From the plot it can be seen that a Torb of roughly 
0.1 can be adopted as a general rule-of-thumb in order to allow 
for an effective constant acceleration appr oximation, in agreem ent 
with the analytical sensitivity calculated in iRansom et al.l ( l2003l) . In 
addition, we point out a tendency for the orbital phases with a sig¬ 
nificant a, to prefer shorter Torb and vice versa for the orbital phases 


® Out of the six Keperian parameters, the longitude of ascending node is 
related to the 3D orientation of the binary orbit, which is not relevant in 
pulsar searching. Hence only the reminding five parameters ai'e considered. 
® We note that the two test data sets employed are taken on different dates 
with different instrumental set-ups, it is thus not appropriate to compai'e 
their detected S/N directly. However, one can attempt a qualitative compar¬ 
ison by normalising the highest S/N of any particular orbital phase to unity. 
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Figure 3. Comparison of the effectiveness in mitigating periodic RFI by applying vaiying extents of RFI cleaning. We show the power spectrum created 
from the time series at DM = 0cm~^pc of an actual survey observation, zooming into the region around 150 Hz where RFI due to hannonics of the 
Australian mains power supply can be seen clearly. In each panel, the original uncleaned spectrum is plotted in red, whereas the spectrum after each RFI 
cleaning procedure is plotted in black. In this particular case, the empirically-generated periodic RFI list removed only 0.01 per cent of the Fourier spectral 
bins (panel (d)), but is still more effective than applying the fixed periodic RFI list which always removed 0.57 per cent of the spectrum (panel (c)). We note 
that subtle differences in the spectra are a result of the change in weighting due to normalisation. 



Figure 4. Detected S/N versus ^orb foi' Ihe two test data sets (test data set 1: 
triangle symbol, test data set 2: circle symbol) at the four selected orbital 
phases. Phases 0.25 and 0.75 represent orbital phases where a is largest. 
Phases 0 and 0.5 correspond to orbital phases where d ^ 0. 

with an d ^ 0. The introduction of eccentricity in the orbital mo¬ 
tion will alter this picture. When the line of sight d is significantly 
non-zero, we expect smaller rorb to perform better (i.e., the peaks 
of these curves shift towards the left). Otherwise, in the less accel¬ 
erated part of the eccentric orbit, slightly larger rorb can result in 
a higher S/N. For a qu antification of the d etectability of eccentric 
binary pulsars see e.g.. lBagchi et al.l (l2013h . 

3.2.2 Acceleration ranges 

Assuming a circular orbit, we derive the theoretical maximum or¬ 
bital acceleration (Umax) for a given Porb by applying Kepler’s 


third law. We consider the upper limit case by setting the orbital 
inclination to be edge on (i.e., i = 90°): 



where Tq = GMq/c^ = 4.925490947 [is, x is the projected 
semi-major axis of the pulsar orbit and c is the speed of light. The 
mass function / is defined by 

(mp + 

where sinz = 1 and the pulsar mass, nip, is taken to be 1.4 Mq. 
The companion mass, rric, then remains the only variable. Fig.|5] 
shows the theoretical maximum orbital acceleration as a function 
of orbital period. We have plotted three scenarios corresponding to 
binary systems with a 0.2 Mq helium white dwarf companion (NS- 
WD), with a second 1.4 Mq neutron star companion (NS-NS) and 
with a hypothetical black hole companion of mass 10 Mq (NS- 
BH). We overplot all published relativistic pulsar binary systems 
with Porb less than 12 hr and significant orbital acceleration reach¬ 
ing above ±1 ms“^ as a reference. 

Based on Fig.[5l sensible acceleration ranges (Aa) can be de¬ 
termined for any particular orbital period. We note that the effect 
of orbital eccentricity has not been taken into account. Flighly- 
eccentric relativistic binary systems can overshoot these theoreti¬ 
cal curves of maximum acceleration significantly, in particular dur¬ 
ing the orbital phase of periastron. Notable examples in Fig. are 
two NS-NS systems PSRs B1913-F16 and B2127-fllC with or¬ 
bital eccentricities of 0.62 and 0.68 respectively. To determine the 
theoretical maximum acceleration while including eccentric sys¬ 
tems would dramatically increase the acceleration search parame¬ 
ter space, making the data processing unfeasible given the current 
computing resources available. Fortunately, a simple consideration 
of the Kepler’s third law shows that, these eccentric systems spend 
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Orbital period (hr) 


Figure 5. Maximum orbital acceleration versus orbital period assuming circular orbits for binary systems of NS-WD (blue line), NS-NS (red line) and NS-BH 
with a 10 Af 0 BH (black line). The maximum orbital acceleration for all published relativistic binaries are also shown as a reference. The coloured regions 
correspond to the parameter spaces prohed by different segments as explained in Section [3.2.3l 


only a relatively brief moment during the highly-accelerated orbital 
phase near periastron, whereas the majority of their orbital phases 
are confined within the less accelerated regime. Therefore, we jus¬ 
tify that our acceleration ranges remain a reasonable compromise 
given the computing resources available. 


3.2.3 Partially-coherent segmentation 

As established in Section [3.2.1l the constant acceleration approxi¬ 
mation is most effective at ro^b ~ 0.1. Hence in order to remain 
sensitive to a wide range of binary orbits it is strategic to use differ¬ 
ent integration lengths. Taking this as a rule-of-thumb, the 4300 s 
full-length observation (s = 1) of the HTRU Galactic plane sur¬ 
vey would enable us to detect binary systems with Porb ^ 12 hr. 
Halving the observation into two equal segments (s = 2) will cor¬ 
respond to binaries of Porb ^ 6 hr, quartered-length observation 
(s = 4) will correspond to Porb ^ 3 hr and segmenting our obser¬ 
vation to one-eighth (s = 8) will correspond to Porb ^ 1-5 hr. We 
set aside binary systems with shorter Porb for future re-processing. 
This is because given the long coherent integration length of this 
survey, such tight binaries should be most efficientl y detected us¬ 
ing th e ‘phase modulation’ technique described in iRansom et al.l 
( l2003h . 

We search for binary systems using an acceleration range 
appropriate for the length of each segment, as discussed in Sec¬ 
tion The maximum acceleration attainable from a IOMq 

NS-BH system is of the order of 1200 m s“^ with a Porb of 1.5 hr, 
hence we have adopted this as the upper bound of our acceleration 
range for the shortest s = 8 segments. The corresponding param¬ 
eter space probed is shown as the pink region in Fig. In order 
to maximise our detectability towards a NS-BH binary, we addi¬ 
tionally search a Aa between ±(200 — 500) ms“^ with the s = 4 
segments (providing sensitivity in the blue region), as well as a Aa 
between ±(1 — 200) ms“^ with the s = 2 segments (providing 
sensitivity in the green region). For the full length s = 1 observa¬ 
tion we search a Aa of ±1 m s“^, which should allow the detection 
of all mildly accelerated binary systems with large orbital periods 
of Porb -S 12 hr, as well as all isolated pulsar systems. 


A schematic of the final pipeline is illustrated in Fig.|6l We call 
it the ‘partially-coherent segmented acceleration search’, as each 
segment is analysed coherently while across the segments the re¬ 
sults are interpreted independently. We search the four configura¬ 
tions (s = 1, 2, 4, 8) in parallel, which is essentially equivalent to 
conducting a multiple pass survey. On one hand, an additional ad¬ 
vantage of this partially-coherent scheme is that by independently 
analysing segments we are less susceptible to some epochs where 
the detection of the pulsar is more difficult, for instance at an or¬ 
bital phase where a is significantly non-zero in a highly-eccentric 
orbit, an intermittent pulsar in switching off phase or a scintilla¬ 
tion induced reduction in S/N. On the other hand, we note that this 
scheme relies on the fact that any binary pulsar to be found has to 
be detectable in at least one of the shortened segments. To be able 
to coherently combine the acceleration search results across seg¬ 
ments would allow us to exploit the full sensitivity achievable with 
the deep integration of this survey, detecting even the weakest rel¬ 
ativistic binary systems. This is the primary goal of our future data 
re-processing. 


4 RE-DETECTIONS OF KNOWN PULSARS 


To verify that no previously known pulsar has been missed by the 
survey, we compute the expected S/N of every pulsar for any par¬ 
ticular observation using the radiometer equation: 


S/N,,,, 


f GVrip/intA/ ^ ^ 


(5) 


The first fraction of Equation © contains parameters related to the 
observational set-up. The ‘degradation factor’ /3 is due to digitisa¬ 
tion and is ~1.16 for our case. The system temperature, Tsys, in¬ 
cludes contribution from both the sky temperature (Tsky ~ 7-6 K) 
and the receiver temperature (Tree ~ 23 K), while G is the antenna 
gain in K Jy“^ and its value varies between 0.58 1, 0.690 and 0.73 5 
depending on the receiver beam (see Table 3 in iKeith et alj|20l3) . 
The number of polarisations summed, rip, is always 2 in our case. 
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Figure 6. A schematic diagram for the ‘partiaiiy-coherent segmented’ pipeiine adopted for this survey. As we progressiveiy downsampie the observation after 
2 X the diagonai DM (DM = i22 cm“® pc) is reached, we end up with six groups of DM ranges with different tsamp (iiiustrated by the coiour density in this 
schematic) each with different number of DM triais (triaisjjM)- Each of the subsequent time series is then passed to the four iower paneis (coiTesponding to 
the four configurations of s = 1, 2,4, 8) to be Fourier transformed, where the number of sampies used in each Fourier transform is marked as ripF"!" We 
show aiso the acceieration range, Aa, the acceieration step size, (5a, as weii as the resuitant number of acceieration triais per DM, triaiSacc/DM. 


The integration time, fint, is 4300 s and / is the effective bandwidth 
of the receiver which is 340 MHz. 

The remaining terms are pulsar dependent and we have taken 
the published data for these terms from PSRCAT. These parameters 
include the spin periods, P, the pulse widths at 50per cent, W 50 , as 
well as the flux densities as observed at 1.4 GHz. Note that we take 
into account the relationship between the reduction in expected flux 
density, Sexp, and the catalogue flux density, S'rroo^ if the pulsar is 
offset from the beam centre. This is given by 

Sexp — >5*1400 ^ J J ! (b) 


where 0 is the radial distance between the published pulsar position 
and the centre of the relevant beam. Assuming a Gaussian drop-off 
of beam sensitivity with a FWHM of 14'.4, we define a to be: 


_ FWHM 
^ ^ 2V2\n2 


(7) 


The mini mum detectable S/N is bas ed on the false alarm statis¬ 
tics (see e.g.. iLorimer & Kramej[200^ and for this survey it is 
taken to be 9. Within the 50 per cent of processed HTRU Galac¬ 
tic plane observations, we record 912 independent known pulsar 
re-detections from 435 pulsars (refer to a re-detection table avail¬ 
able as online material, a sample of which can be see in Table 1C ll 
in the Appendix). 

Fig.|7]plots the observed S/Nobs versus the expected (S/Nexp) 
calculated from Equation © for the known pulsar re-detections. 
Out of the 912 re-detections, 198 lie outside the 14\4 FWHM of the 



Figure 7. A comparison of the observed S/Ns of known pulsai' re-detections 
with their predicted values. The solid line shows the expected 1:1 correla¬ 
tion. The dotted and the dashed lines mark the region containing 95 and 99 
per cent of the data points. Pulsars with DM^lOOcm”^ pc are potentially 
affected by scintillation and are coloured in green. The eiTor bars reflect the 
published uncertainties in the pulsai* flux densities from PSRCAT. 
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receiver beam. Note that the sensitivity of the telescope outside the 
FWHM of the primary beam is complicated hy its sidelohe pattern 
and is not well quantified. Hence these 198 re-detections are disre¬ 
garded for the purpose of S/N comparison. Another five detections 
come from known pulsars with no published catalogue flux den¬ 
sity. They are therefore not included in Fig. |7] The black points in 
Fig.|7]denote the remaining 709 known pulsar re-detections within 
the FWHM. It can he seen that the correlation between S/Nexp and 
S/Nobs is very good. Most of the data points lie close to the line of 
1:1, implying the wider A/ and the longer tint of the HTRU survey 
has provided the increa sed sensitivity expected. 

iKeith et ^ (1201 Oh have pointed out that data points below the 
1:1 correlation could he influenced by the bias resulting from the 
tendency to publish the discovery observation which is more likely 
to have the highest S/N due to scintillation, leading to the higher 
catalogue flux densities. Despite this, a few re-detections are sig¬ 
nificantly worse than expected and lie below the 99 per cent lower 
boundary. They are likely to be due to scintillation and/or contam¬ 
ination by RFI. In the case of the very bright pulsars on the top 
right-hand comer of Fig. [T] the loss in survey sensitivity might be a 
consequence of the 2-bit digitisation. This, however, is not a prob¬ 
lem for our purpose here. 

Out of all the expected known pulsar re-detections, 22 have 
been missed in the data processing (refer to Table IC2I in the Ap¬ 
pendix for potential reasons of these non-detections). Folding the 
observation directly with the catalogue pulsar ephemerides resulted 
in the recovery of seven pulsars (listed as S/Neph in Table IC2l l. 
Among the missed re-detections, one is attributed to the magne- 
tar PSR 11622—4950. This magnetar is being followed up by an 
on-going HTRU timing programme. However, its radio emission 
appears to have ceased since October 2010, which is consistent 
with the non-detection of the HTRU Galactic plane survey obser¬ 
vation taken in December 2010. The radio emission of magnetars 
are known to fade with time, hence we have likely observed the 
end of the radio emission of PSR 11622—4950. Excluding several 
other similar cases (refer to Appendix [C^ for the reasoning) where 
phenomena intrinsic to the pulsar have prevented the re-detections, 
we conclude that only ~ 1 per cent of the published pulsars ex¬ 
pected have been missed and the large number of re-detections in¬ 
dicate that the HTRU Galactic plane survey is performing as ex¬ 
pected with little loss of sensitivity. We also remark that none of 
the missed pulsars belongs to binary systems (refer to Table lC3l in 
the Appendix). 


5 NEWLY-DISCOVERED PULSARS 

A total of 60 pulsars have been discovered in the 50 per cent pro¬ 
cessed data of the Galactic Plane survey. Listed in Table [2] are 
their folded signal-to-noise ratios (S/Nhtru) at discovery. The 
previous PMPS as well as the medium-latitude sub-survey of the 
HTRU have a complete overlap in the region of sky with the HTRU 
Galactic plane survey. We have inspected the PMPS and medium- 
latitude archival data to determine if any of the newly-discovered 
pulsars were detectable and listed the respective folded S/Npmps 
and S/Nmediat in columns 3 and 4 of Table|^ Eighteen and five pul¬ 
sars have respective S/Npmps and S/Nmediat greater than the the¬ 
oretical detection threshold of S/Nmin.PMPS =S/Ninin,medlat = 8, 
which means they could have been discovered in these archival ob¬ 
servations. They might have been missed due to the large number of 
candidates produced in the PMPS and medium-latitude processing. 
Improvements in the RFI mitigation scheme of the current HTRU 


Table 2. Table listing the S/N, flux density (51400), derived luminosity 
(Li4oo) and pulse widths (IV50, Wio) of the 60 newly-discovered pulsars. 


PSR Jname 

S/N 

HTRU 

S/N 

medlat 

S/N 

PMPS 

5i400 

(mJy) 

W 50 

(ms) 

IFio 

(ms) 

-^1400 
(mJy kpc^) 

J1002-5919 

18.2 

6.8 

7.7 

0.17 

61.3 

69.1 

12.6 

J1101-6424 

19.7 

7.4 

- 

0.27 

0.54 

1.39 

5.4 

J1151-6108 

10.0 

- 

- 

0.06 

5.65 

9.89 

1.1 

J1227-63 

19.4 

— 

— 

— 

15.0 

34.4 

— 

J1244-6359 

12.5 

- 

- 

0.15 

9.8 

22.9 

A.l 

J1248-6444 

13.4 

— 

< 6.4 

0.15 

54.3 

78.5 

3.3 

J1255-62 

14.2 

- 

- 

0.13 

34.8 

65.0 

21.8 

J1349-63 

17.4 

— 

9.7 

— 

23.0 

49.7 

— 

J1525-5523 

12.6 

10.3 

10.5 

0.21 

31.0 

40.9 

1.1 

J1528-5547 

14.4 

- 

- 

0.07 

52.9 

110 

1.2 

J1532-56 

12.0 

— 

9.7 

0.10 

27.0 

54.0 

1.8 

J1538-5621 

22.7 

9.2 

15.7 

0.14 

34.1 

60.6 

1.2 

J1612-49 

12.6 

- 

- 

0.16 

28.7 

73.1 

10.3 

J1612-55 

17.7 

— 

— 

0.11 

37.6 

56.5 

3.6 

J1616-5017 

21.5 

— 

— 

0.17 

9.76 

20.5 

3.8 

J1622-4845 

12.1 

- 

- 

0.15 

17.5 

42.3 

3.1 

J1627-49 

14.1 

— 

8.6 

0.13 

69.3 

97.1 

6.0 

J1627-51 

21.6 

- 

13.3 

0.08 

9.8 

16.3 

0.98 

J1634-49 

14.3 

— 

8.5 

0.13 

43.5 

62.5 

10.5 

J1638-44 

15.5 

- 

- 

0.15 

31.6 

50.6 

6.8 

J1649-3935 

17.2 

6.8 

9.6 

0.05 

22.8 

45.7 

1.6 

J1658-47 

11.5 

6.9 

7.7 

0.18 

27.4 

43.8 

23.0 

J1708-3641 

13.5 

- 

< 6.5 

0.12 

21.8 

78.3 

2.7 

J1710-37 

13.3 

— 

< 7.7 

0.10 

61.4 

105 

2.2 

J1718-41 

12.4 

- 

- 

0.07 

28.4 

69.1 

2.0 

J1720-36 

14.5 

— 

— 

0.08 

8.41 

10.8 

1.4 

J1723-38 

14.3 

- 

10.6 

0.08 

13.7 

35.0 

2.8 

J1730-34 

15.4 

- 

7.1 

0.03 

25.9 

51.7 

1.5 

J1731-3322 

10.5 

— 

< 6 

0.10 

33.5 

71.4 

12.5 

J1732-35 

11.1 

- 

- 

0.13 

6.6 

14.1 

3.3 

J1734-3058 

21.0 

- 

9.3 

0.11 

10.7 

22.6 

1.6 

J1738-2736 

26.0 

— 

9.5 

0.17 

13.7 

23.7 

4.4 

J1741-34 

12.0 

— 

— 

0.20 

19.1 

43.4 

3.4 

J1743-35 

14.4 

— 

9.7 

0.05 

16.9 

33.8 

0.5 

J1746-27 

8.7 

- 

- 

0.15 

29.9 

42.6 

5.6 

J1748-30 

12.3 

6.0 

8.9 

0.16 

25.5 

79.4 

10.5 

J1750-28 

11.6 

- 

- 

0.09 

16.8 

35.1 

2.3 

J1755-25 

15.1 

— 

9.0 

0.14 

11.7 

25.7 

14.9 

J1755-26 

15.6 

6.9 

— 

0.14 

16.3 

25.7 

3.8 

J1756-25 

20.5 

- 

- 

0.20 

25.4 

38.0 

18.4 

J1757-27 

14.6 

— 

— 

0.07 

0.49 

0.88 

2.0 

J1759-24 

28.5 

7.6 

- 

0.50 

72.3 

159 

59.4 

J1811-1717 

14.2 

6.4 

< 5.8 

0.20 

38.8 

59.0 

9.0 

J1819-17 

11.2 

— 

— 

0.14 

87.1 

157 

0.5 

J1824-1350 

14.6 

- 

- 

0.08 

19.4 

58.2 

6.5 

J1825-1108 

17.5 

— 

7.5 

0.13 

42.1 

68.8 

0.9 

J1829-1011 

19.8 

- 

9.3 

0.25 

98.8 

125 

11.2 

J1830-10 

9.9 

— 

< 7.3 

0.10 

5.36 

12.7 

1.3 

J1833-0209 

11.9 

— 

< 6.8 

0.06 

17.4 

30.7 

3.1 

J1834-09 

18.1 

- 

- 

0.17 

64.5 

152 

7.0 

J1835-0924b 

13.0 

8.9 

— 

— 

27.1 

45.3 

— 

J1835-0928 

10.5 

- 

- 

- 

39.5 

64.2 

- 

J1836-11 

10.8 

- 

- 

0.12 

25.7 

55.1 

3.2 

J1838-0107 

27.1 

10.5 

9.7 

0.05 

34.4 

45.0 

1.9 

J1839-0223 

11.2 

- 

8.3 

0.11 

45.3 

67.9 

4.5 

J1839-0332 

20.2 

8.6 

12.7 

- 

31.9 

63.7 

- 

J1842-0800 

9.8 

— 

9.9 

0.10 

24.9 

52.3 

1.7 

J1843-0510 

11.8 

- 

7.5 

0.07 

24.0 

36.0 

1.9 

J1844-0302 

13.6 

— 

— 

0.12 

19.0 

33.3 

6.4 

J1847-0427 

11.4 

- 

- 

- 

30.9 

48.4 

- 
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Table 3. Specifications of the observing system employed for the timing 
observations in this work. The antenna gain is represented by G and Tsys 
is the receiver system temperature. The central frequency in MHz is repre¬ 
sented by /c and B is the bandwidth in MHz. 


Receiver 

G 

(KJy-i) 

Tsys 

(K) 

Backend 

/c 

(MHz) 

B 

(MHz) 

Multibeam 

0.74 

23 

Parkes DFBs 

1369 

256 




Parkes BPSR 

1352 

340 




Parkes APSR 

1369 

256 

Single-pixel 

1.00 

28 

Lovell DFB 

1532 

384 




Jodrell ROACH 

1532 

400 


processing as mentioned in Section[3T|have likely helped avoiding 
such candidate confusion, enabling the detection of these pulsars. 
A further twelve and eight pulsars have weak detections from these 
archival observations, with respective S/Npmps and S/Nmediat less 
than 8. In the cases when only tentative detections are suggestive 
we have listed an upper limit for their S/Npmps or S/Nmediat in 
Table|2 These pulsars might have been missed in the processing of 
these archival surveys as their low S/Ns might have prevented them 
from being selected for visual inspection. The remaining 25 pul¬ 
sars are detectable from neither the PMPS nor the medium-latitude 
archival data. 

Full flux density calibration is implemented for 54 pulsars 
with observations taken with the Parkes Digital Filter bank sys¬ 
tems (DFB). We calibrate each observation by using an averaged 
observation of Hydra A, and we account for the differential gain 
and phase between the feed with an observation of the noise diode 
coupled to the the feeds. It is important that this calibration is taken 
adjacent to the targeted pulsar observations. Column 5 of Table 
reports the mean flux densities averaging over all available timing 
observations for each pulsars at 1.4 GHz, Si 4 oo- In turn, we infer 
the luminosity at 1.4 GHz, L 1400 = <S'i 4 oo x cP, where d is the dis¬ 
tance of the pulsar in kpc (see column 9 of Table in accord ance 
to the NE2001 electron density model jCordes & Lazidl2002h . 

Also listed in columns 6 and 7 of Tablej^are the pulse widths 
measured at 50 (IT 50 ) and 10 (PFio) per cent of the highest peak. 
The pulse profile of each pulsar at 1.4 GHz is shown in Fig. [8] For 
the pulsars with coherent timing solutions, we summed all obser¬ 
vations to form high S/N mean profiles. Otherwise, we plot the 
profiles from the single observation at discovery. The peak of the 
profiles have been normalised to unity and their peaks placed at 
phase 0.2. Mo st of the long-period normal pulsars have typical 
pulse profiles jLvne & Graham-Smithll200^ . with single-peaked 
pulses and 5 < 10 per cent. One of the 60 newly-discovered pul¬ 
sars, PSR J1847—0427, shows a broad pulse profile with a hint of 
an interpulse and is further discussed in S ectionl6.6l These results 
are en tirely consistent with the findings of IWeltevrede & Johnstorj 
( l2008h . that ~2 per cent of the known pulsar population is observed 
with an interpulse. A few of the pulsars with high DM display the 
classical exponential tail of scattering caused by propagation of the 
radio signal through the interstellar medium. 

When a pulsar is first discovered, our knowledge of its sky 
position, rotation period and DM are only approximate. Follow up 
timing observations of at least 1 yr are necessary to precisely de¬ 
termine its rotational, astrometric and, if any, orbital parameters. 
We sum each timing observation over both frequency and time to 
produce an integrated pulse profile. The PSRCHIVE data analysis 
package jHotan et alj|2004n is used to convolve a noise-free ana¬ 
lytic reference templa te with each individual profile to produce a 
time of arrival (TOA; iTavlodfl?^ . The TEMP02 software pack¬ 


age dHobbs et alj|200^ is then employed to fit a timing model 
to all TOAs of the pulsar. Towards the end of the timing analy¬ 
sis procedure when the respective reduced comes close to one, 
we can then assume a reliable fit is achieved which is only influ¬ 
enced by the presence of radiometer noise in the template. As a 
last step, we compensate for these systematic effects by calculating 
dataset-specific calibration coefficients (also known as ‘EFAC’ in 
TEMP02). These coefficients are applied to scale the TOA uncer¬ 
tainties such that each final respective reduced is unity. 

All observations presented here have been taken at 1.4 GHz at 
Parkes using backends including the DEBs with incoherent dedis¬ 
persion and the ATNF Parkes Swinburne Recor defl (APSR) with 
coherent dedispersion dHankins & Ricke3l975h . Pulsars with dec¬ 
lination North of —35^ were timed at the Jodrell Bank Observa¬ 
tory with the Lovell 76-m telescope, using a DFB backend and a 
ROACH backend. The latter is based on the ROACH FPGA pro¬ 
cessing boarcfl and coherently dedisperses the data. Refer to Ta- 
blej^for the specifications of all observing systems employed. Tim¬ 
ing observations at Jodrell Bank were performed approximately 
once every three weeks, whereas Parkes observations are more ir¬ 
regular with gaps ranging from days to months depending on tele¬ 
scope availability. At both telescopes, integration times of at least 
20 min are typically required, with longer observations for weaker 
pulsars to achieve adequate S/N of at least 10. 

Table|4]presents the 30 newly-discovered isolated pulsars with 
coherent timing solutions, whereas the three newly-discovered bi¬ 
nary systems are listed in Table [5] The Damour-Deruelle (DD) 
timing model dPamour & Deruellari986h in TEMP02 is a theory- 
independent description for eccentric binary orbits. However, for 
binaries with small eccentricities the location of periastron is not 
well-defined and using the DD timing model results in a high co- 
variance between the longitude of periastron (at) and the epoch 
of peri astron (Tp). For thse pulsars, we use the ELLl timing 
model jLange et al.ll200ll) alternatively. The ELLl timing model 
avoids the covariance by using the Laplace-Lagrange parameters 
(ei = e sin uj and 62 = e cos at) and the time of ascending node 
passage (Tasc) instead of To as in the DD timing model. A re¬ 
maining 28 of the most recent discoveries have not yet been allo¬ 
cated sufficient follow-up telescope time, but will be monitored in 
the coming months. For these pulsars we have reported their dis¬ 
covery parameters in Table [6] 

No associated gamma-ray pulsations have been identified for 
any of the newly-discovered pulsars presented in this paper. Cur¬ 
rently, the short time span of the radio ephemerides are inadequate 
to phase-fold the now more than five years data from the Large Area 
Telescope (LAT) on the Fermi Gamma-Ray Space Telescope. Con¬ 
tinuous follow-up timing of these pulsars and any future discoveries 
from the HTRU Galactic plane survey would provide up-to-date ra¬ 
dio ephemerides critical for recovering any associated gamma-ray 
pulsations. We note that for the 32 pulsars currently with P and P 
meas urements, none of them has high log vT/(see lAbdo et Tl 
I 2 OI 3 L for details), hence unlikely to be detected by Fermi in the fu¬ 
ture. 


^ http://astronomy.swin.edu.au/pulsar/?topic=apsr 

® In practice, for the declination between —30° and —35°, only bright 
pulsars ai'e followed-up at the Lovell Telescope. This is because of the short 
visible hours, as well as the challenges associated with such low elevation 
observations, namely the need for low wind conditions and the potential 
high spillover. 

® https;//casper.berkeley.edu/wiki/ROACH 
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Pulse Phase 

Figure 8. Average pulse profile of the 60 newly-discovered pulsars, plotted by placing the peak at phase 0.2 and scrunching to 128 phase bins. Also marked 
on each profile is the respective pulsai* names, the spin periods in s followed by the DMs in cm~^ pc. 





































































Table 4. TEMP02 best-fitting pai'ameters of 30 newly-discovered isolated pulsars from the HTRU Galactic plane survey. We list the equatorial (R.A. and Dec.) and Galactic (Z and b) position, the spin period and the 
DM of these pulsars. Values in parentheses are the nom inal 1 a uncertain t ies in the last digits. Pulsars for which no full timing solution is available have been assigned a temporary name containing only two digits 
of declination. Pulsar distances are derived according to lCordes & Lazicl i2Q02h . We include fitting related parameters such as the data span, the reference epoch, the number of TOAs employed (Ntoa the RMS 
of the TEMP02 fit and the reduced ,). The characteristic age (Tr-). the surface magnetic field and the spin-down energv (E) are derived using equations which can be found inlLorimer & Krameil 

m . 


K) 


PSR name 

R.A. (J2000) 

^hms^ 

Dec. (J2000) 

r'") 

1 

n 

b 

n 

P 

(ms) 

P 

(10-1®) 

DM Dist Data span 

(cm-®pc) (kpc) (MJD) 

Epoch 

(MJD) 

N 

TOA 

RMS 

(its) 

Xred* 

Tc 

(Myr) 

^surf E 

(lO^^G) {lO^^evg 

J1002-5919 

10:02:20.50(5) 

-59:19:37.2(6) 

282.83 

-3.22 

713.4892535(6) 

104(10) 

347(2) 

8.7 

56499-56930 

56677 

28 

1262 

1.3 

110 

28 

14 

J1151-6108 

11:51:56.86(2) 

-61:08:17.6(3) 

295.81 

0.91 

101.633196293(2) 

10278.4(10) 216(1) 

4.2 

56499-56930 

56714 

23 

696 

1.5 

0.16 

100 

390000 

J1248-6444 

12:48:32.87(2) 

-64:44:00.0(8) 

302.62 

-1.86 

1234.89349389(3) 

1932(16) 

236.4(11) 

4.8 

56498-56931 

56715 

14 

567 

0.6 

10 

160 

41 

J1525-5523 

15:25:36.064(7) 

-55:23:27.2(2) 

323.66 

1.15 

355.1560369998(16) 

8.96(7) 

124.7(3) 

2.3 

54792-56571 

55682 

29 

396 

0.6 

630 

5.7 

7.9 

J1528-5547 

15:28:39.18(12) 

-55:47:23(2) 

23.80 

0.58 

3467.3015170(17) 

7750(20) 

269 

4.2 

56127-56967 

55839 

27 

7798 

8.2 

7.1 

525 

7.4 

J1538-5621 

15:38:43.29(2) 

-56:21:55.5(4) 

324.61 

-0.70 

1908.49449143(2) 

7311.8(7) 

150(18) 

3.0 

54939-56886 

55912 

14 

926 

0.9 

4.1 

380 

42 

J1616-5017 

16:16:29.89(4) 

-50:17:14.9(2) 

332.83 

0.29 

491.384121830(13) 

46592(3) 

194(4) 

4.7 

56469-56933 

56701 

14 

597 

3.4 

0.17 

480 

16000 

J1622-4845 

16:22:05.45(10) 

-48:46:29(4) 

334.53 

0.73 

735.09164794(8) 

110 (20) 

292(7) 

4.6 

56499-56968 

56715 

12 

1056 

1.3 

98 

30 

12 

J1649-3935 

16:49:06.68(3) 

-39:35:44.2(14) 

344.58 

3.34 

770.909760511(14) 

39.1(4) 

290(7) 

5.4 

55085-56886 

55985 

15 

1195 

1.5 

310 

18 

3.4 

J1708-3641 

17:08:35.91(3) 

-36:41:21.5(14) 

349.23 

2.13 

587.566793839(14) 

134.9(5) 

316(3) 

4.8 

55208-56886 

56047 

19 

1581 

0.6 

69 

29 

26 

J1731-3322 

17:31:14.31(10) 

-33:22:45(2) 

354.60 

0.22 

544.67054719(4) 

28093(18) 

877.4 

11.0 

56500-56932 

56716 

26 

2931 

1.5 

0.3 

400 

6900 

J1734-3058 

17:34:50.839(7) 

-30:58:41.4(8) 

357.03 

0.89 

541.285680730(5) 

16.2(5) 

260.0(14) 

3.8 

55976-56784 

56400 

128 

1750 

2.8 

530 

9.5 

4.0 

J1738-2736 

17:38:14.628(5) 

-27:36:25.8(9) 

0.27 

2.08 

627.715518484(19) 

4624.3(7) 

323.6(7) 

5.1 

56322-56884 

56315 

69 

760 

3.6 

2.2 

170 

740 

J1750-28 

17:50:04.0(2) 

-28:45(7) 

0.66 

-0.74 

1300.5131462(7) 

5780(190) 

388(12) 

5.1 

56681-56884 

56783 

22 

1567 

1.7 

3.6 

280 

100 

J1755-26 

17:55:16.29(4) 

-26:00(1) 

3.62 

-0.33 

430.87215605(14) 

12280(50) 

405(4) 

5.2 

56684-56940 

56827 

30 

1207 

2.4 

0.55 

230 

6100 

J1757-27 

17:57:54.795(2) 

-27:45(7) 

2.40 

-1.72 

17.68721477177(11) 

0 .21(2) 

334 

5.3 

56675-56881 

56778 

50 

24 

2.0 

1300 

0.20 

1600 

J1811-1717 

18:11:26.53(19) 

-17:17:47(14) 

13.04 

0.69 

391.38510099(9) 

120(80) 

545.5 

6.7 

56640-56884 

56762 

39 

2461 

1.3 

49 

22 

83 

J1824-1350 

18:24:50.181(17) 

-13:50:21.0(16) 

17.61 

-0.52 

1396.59854604(4) 

607(7) 

551(7) 

6.5 

56341-56884 

56612 

75 

2957 

1.9 

36 

93 

8.8 

J1825-1108 

18:25:18.3(2) 

-11:08:54(12) 

20.05 

0.64 

1925.8717727(8) 

1900(600) 

121(19) 

2.7 

56682-56884 

56783 

29 

1003 

3.5 

16 

200 

11 

J1829-1011 

18:29:06.06(4) 

-10:11:25(4) 

21.33 

0.26 

829.16602078(16) 

-5(39) 

610 

6.7 

56626-56980 

56755 

44 

2962 

0.7 

> 170+ 

< 25+ 

< 5.2+ 

J1833-0209 

18:33:05.411(8) 

-02:09:16.4(3) 

28.92 

3.09 

291.930632181(4) 

2754.8(7) 

325.4 

7.2 

56341-56886 

56614 

97 

1535 

1.8 

1.7 

91 

4400 

J1835-0924b* 

18:35:21.822(16) 

-09:24:15.9(10) 

22.74 

-0.74 

235.248924972(4) 

11.84(3) 

500 

6.3 

55905-56885 

56394 

142 

4576 

2.7 

315 

5.3 

36 

J1835-0928 

18:35:22.22(3) 

-09:28:02(2) 

22.68 

-0.77 

621.73399423(3) 

969(2) 

450 

5.8 

55906-56885 

56395 

98 

6772 

5.9 

10 

79 

160 

J1838-0107 

18:38:39.423(2) 

-01:07:48.64(10) 

30.48 

2.33 

444.4257246169(17) 

5.46(15) 

268.9 

6.1 

56094-56884 

56489 

131 

602 

2.2 

1300 

5.0 

2.5 

J1839-0223 

18:39:58.032(3) 

-02:23:09.5(16) 

29.50 

1.46 

1266.79012424(5) 

4763(6) 

330 

6.4 

56094-56881 

56488 

117 

9103 

6.3 

4.2 

250 

92 

J1839-0332 

18:39:56.583(12) 

-03:32:58.6(6) 

28.46 

0.93 

2675.68226451(6) 

4760(6) 

195.1 

4.8 

56098-56885 

56491 

101 

2820 

1.7 

8.9 

360 

9.8 

J1842-0800 

18:42:54.943(13) 

-08:00:53.5(7) 

24.83 

-1.77 

1255.46857429(3) 

194(5) 

188.6 

4.1 

56341-56880 

56611 

57 

1949 

0.9 

100 

50 

3.8 

J1843-0510 

18:43:09.732(12) 

-05:10:04.5(5) 

27.39 

-0.52 

671.613828819(13) 

3892.3(12) 

257 

5.2 

56144-56885 

56514 

98 

2572 

2.0 

2.74 

160 

510 

J1844-0302 

18:44:06.918(12) 

-03:02:11.2(5) 

29.40 

0.24 

1198.63027492(2) 

7809(5) 

533 

7.3 

56341-56885 

56613 

87 

2586 

1.3 

2.4 

310 

180 

11847-0427 

18:47:18.86(5) 

-04:27:59(2) 

28.49 

-1.12 

259.24664288(2) 

5.9(5) 

188.3 

4.7 

56341-56882 

56611 

71 

9527 

3.3 

69 

13 

140 


* Note that PSR J1835—0924b is unrelated to the previously known PSR J! 835—0924, which has a catelogue spin period of 859.192 ms at a DM of 471 cm~^ pc. 

^ For PSR J1829—1011, the period derivative-related parameters are derived with the 2a upper limit of P < 7.4 x 10“^^. 

^ The reduced stated here represents the value before the application of EFAC. Note that the rest of the timing solutions have EFACs incorporated, bringing the reduced fo unity. 
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Figure 9. Spin period versus minimum companion mass for all published 
pulsars in binaiy sys tems. We have clas sified their companions according 
to the description in lXauris'^ iDiMi- Green points indicate ultra-light 
(UL) binaries, red points indicate helium white dwarf (He-WD) compan¬ 
ions, blue points indicate CO-WD, puiple points indicate main-sequence 
star (MS) companion and grey points indicate NS-NS systems. The three 
newly-discovered binaries are additionally represented by star symbol, and 
only for PSR J1755—25 the companion type is yet undetermined. 


onto the neutron star. Furthermore, given the similarities between 
the two systems (see Table O, we can deduce that the IMXB 
progenitor system of PSR J1101—6424 is likely to have a sim¬ 
ilar initial donor mass of ^ 4.5 dz 0.5 Mm, as e stimated for 
PSR J1614 -2230 jTauris et al.ll201lULin et Sll201 ih . According 
to Fig. 9 of iTauris et alj ( l20llF . this would suggest a neutron star 
mass >1.7 Mq for PSR J1101—6424. Any potential detection or 
constraint of Shapiro delay in PSR J1101—6424 thus implies good 
prospects for measuring a pulsar with high mass, and continued 
monitoring of this binary system is of great interest. 


6.2 PSR J1244—6359, a mildly-recycled binary in a 17.2-d 
orbit 

This system is in a binary with orbital period of 17.2 d and has 
a companion with a median mass of 0.69 Mq (Table O. The 
small orbital eccentricity of 0.000179(19) makes it unlikely that 
PSR J1244-6359 is a DNS system. Rather, PSR J1244-6359 
is most probably a mildly-recycled binary system with a heavy 
CO-WD companion. The spin period of 147.3 ms is on the slow 
side. However, Fig. [9] shows it is not an outlier when compar¬ 
ing to the whole CO-WD binary population. The location of 
PSR J1244—6359 on a P-P diagram is intriguing, bridging the 
population of normal and recycled pulsars (see Fig. d. 


Despite the improved acceleration search algorithm used in 
this analysis, no previously unknown relativistic binary pulsar has 
yet been found. However, all known relativistic binary pulsars in 
the survey region of the 50 per cent processed observations have 
been re-detected with a higher significance than obtained in previ¬ 
ous analyses (refer to Table lC3l in the Appendix for details). 


6 INDIVIDUAL PULSARS OF INTEREST 

6.1 PSR JIIOI—6424, a Case-A Rocbe-lobe overflow cousin 
ofPSRJI6I4-2230 

PSR JllOl—6424 is a binary MSP with an orbital period of 9.6d 
and a spin period of 5.1 ms (Table [Sj. Assuming an orbital incli¬ 
nation of i = 90° and a pulsar mass of 1.35 Mq, we find the 
minimum companion mass of PSR JllOl—6424 to be 0.47 Mq. 
A heavy companion mass as such tends to point to an evo lution¬ 
ary tr ack from an intermediate-mass X-ray binary flMXB: rTaurisl 
I 2 OI ih . with the companion being a carbon-oxygen white dwarf 
(CO-WD). However, the fast spin period of about 5 ms indicates 
a rare example of full-recycling from a long mass-transfer phase. 
This can only be achieved v i a Cas e A Roche lobe overflow (RLO) 
as discussed in iTauris et al.l 1 I 2 OI ih . making PSR JllOl—6424 the 
second known system to descend from this IMXB evolutionary 
track. The only other binary MSP formed from Case A RLO is 
an IMXB the 3.2-ms binary PSR J1614—2230. Fig.|^plots the spin 
periods versus the minimum companion masses of all known pul¬ 
sars as listed in PSRCAT. It can be seen that PSR JllOl—6424 is 
indeed located at the edge of the CO-WD population, together with 
PSR J1614-2230. 

We recall that PSR J1614—2230 is o ne of the heaviest neu - 
tron stars known, with a mass of 2.0 Mq JPemorest et al.ll2010h . 
This does not yet necessarily imply PSR J1101—6424 should also 
have such a heavy pulsar mass. Nevertheless, the long-lasting ac¬ 
cretion phase during Case A RLO leads to significant accretion 


6.3 PSR J1755—25, an eccentric binary pulsar with a heavy 
companion 

PSR J1755—25 was discovered with an unremarkable spin period 
of 315 ms at a DM of 751 cm“® pc. We have now a coherent tim¬ 
ing solution spanning over ~40d and PSR J1755—25 appears to 
be in a binary with an orbital period of 9.7 d and an unusually high 
orbital eccentricity of 0.09 with a heavy companion of minimum 
mass me,min = 0.40 Mq (see Table|5j. Due to the short span of 
the timing data we cannot yet fit for the pulsar position and the pe¬ 
riod derivative. The covariance between these two parameters will 
be resolved when more timing data become available. The nature 
and the binary evolution scenario of this system is interesting and 
will be presented in a future publication. 

6.4 PSR J1757—27, an isolated recycled pulsar 

We have a coherent timing solution for PSR J1757—27 across a 
time span of 557 d with observations taken using the Lovell tele¬ 
scope. PSR J1757—27 has a spin period of 17 ms and a P of the 
order of 10“^®. The small P implies that PSR J1757—27 is not a 
young pulsar. Rather, PSR J1757—27 appears to be an isolated re¬ 
cycled pulsar, or in a very wide binary orbit with P^ih of the order 
of years. Further timing observations is crucial for revealing the na¬ 
ture of this pulsar. If it is indeed proved to be isolated, it will add 
to the currently small population of 30 Galactic (i.e., not associated 
with a Globular Cluster) millisecond-period isolated pulsars. This 
system would have interesting implications to the evolution sce¬ 
narios of millisecond-period isolated pulsars. Currently, the most 
adopted formation scenario of isolated MSPs is that a pulsar sys¬ 
tem with a very-low mass companion cou ld continue ablating it s 
companion until it ceases to exist (see e.g.. lRuderman etn]|l989h . 
However, this would imply fully recycling and hence the pulsar 
should have very fast spin period of the order of a few ms, at odds 
with the 17 ms spin of PSR J1757—27. Another possible scenario 
is that PSR J1757—27 comes from a disrupted binary system of 
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Table 5. TEMP02 best-fit paramete rs for the three newly- discovered binary systems. PSR JllOl—6424 can be compared to that of PSR J1614—2230 listed 
at the last column (data taken from IPemorest et^ j201(t) '). Values in parentheses are the nomin al 1 -cr uncertainties in the last digits. The last panel shows 
derived parameters, th e respective equations for which can be found in IPorimer & Kram^ feOOSh . except for the DM distance which is derived according to 
ICordes & Lazid j2002h . 


Parameter 

11755-25 

11244-6359 

11101-6424 

11614-2230 tDemoresr et al. 20101! 

Right ascension, a (J2000) 

17:55.6(5) 

12:44:47.693(18) 

11:01:37.1923(5) 

16:14:36.5051(5) 

Declination, (5 (J2000) 

-25:53(7) 

-63:59:47.4(3) 

-64:24:39.332(2) 

-22:30:31.081(7) 

Galactic longitude, / (°) 

3.76 

302.20 

291.42 

352.64 

Galactic latitude, b {°) 

-0.34 

-1.13 

-4.02 

20.19 

Spin period, P (ms) 

315.19598238(5) 

147.27431048(2) 

5.109272904279(2) 

3.1508076534271(6) 

Period derivative, P (10“^®) 

_ C 

4.4(8) 

0.0018(2) 

0.0096216(9) 

Dispersion measure, DM (cm~^ pc) 

751(3) 

286(1) 

207 

34.4865 

Orbital period, Porb (d) 

9.69633(2) 

17.170748(2) 

9.6117082(3) 

8.6866194196(2) 

Projected semi-major axis, x (It-s) 

12.2838(2) 

24.0329(2) 

14.02466(3) 

11.2911975(2) 

Epoch of ascending node, Tasc (MJD) 

56904.1259(6)'^ 

56513.0(3)'* 

55689.00791(3) 

52331.1701098(3) 

esinw, ei (10“®) 

_ d 

_ d 

25.6(13) 

0.11(3) 

ecosLJ, 62 (10~®) 

_ d 

_ d 

2.0(12) 

-1.29(3) 

InfeiTed eccentricity, e (10“®) 

0.08932(3)'* 

0.000179(19)'* 

26.0(13) 

1.30(4) 

Longitude of periastron, lj (°) 

129.65(2) 

144(7) 

85(2) 

- 

Minimum companion mass*^, me,min (-^ 0 ) 

0.40 

0.57 

0.47 

- 

Median companion mass^, mc,med (-^ 0 ) 

0.48 

0.68 

0.57 

0.500(6)*' 

Binary model 

DD 

DD 

ELLl 

ELLl 

First TOA (MID) 

56901 

56503 

56401 

52469 

Last TOA (MJD) 

56942 

56968 

56967 

55330 

Timing epoch (MJD) 

56922 

56222 

56597 

53600 

Points in fit 

29 

25 

42 

2206 

Weighted RMS residuals (p.s) 

705 

451 

22 

1.1 

Reduced ^ 

1.4 

1.8 

0.7 

- 

DM distance (kpc) 

10.3 

5.6 

4.5 

1.2 

Characteristic age, Tc (Myr) 

_ C 

520 

44000 

5200 

Spin down energy loss rate, E (10^® ergs“^) 

_ C 

56 

540 

12000 

Chai'acteristic dipole surface magnetic field 
strength, (10^“ G) 

_c 

2.6 

0.0098 

0.018 


“ rric min is calculated for an orbital inclination of i = 90° and an assumed pulsar mass of 1.35 Mq. 

** nic.med is Calculated for an orbital inclination of i = 60° and an assumed pulsar mass of 1.35 Mq. 

° For PSR 11755—25, the time .span of the data is only ~ 0.1 yr, which is not enough to resolve the degeneracy between the position and the spin down rate. 
Hence no period derivate-related par ameters can be derived. 

For PSRs J1755—25 and J1244—6359 the DD model is used. We quote Tq instead of Tasc- e is directly fitted for and not inferred from the e parameters. 

° The reduced stated here represents the value before the application of EFAC. Note that the rest of the timing solutions have EFACs incorporated, bringing 
the reduced x^ to unity. 

^ This companion mass is measured from the detected Shapiro delay. 


Double Neutron Star (DNS). Given that PSR J1757—27 appears to 
have been recycled, it would have to be the first formed pulsar of 
the DNS, which has subsequently been separated from its compan¬ 
ion during the secondary supernova explosion. 

PSR J1757—27 has a small duty cycle of 2.8 per cent. If timing 
residuals with an RMS of the order of a few ps can be achieved 
with future timing observations, PSR J1757—27 can potentially be 
a good candidate for employment in a pulsar timing array. 


6.5 PSR J1759—24, an eclipsing binary system with a long 
orbital period 

PSR 11759—24 was first discovered from a survey observation 
taken at MID 55675.9 with strong, unambiguous pulsar character¬ 
istics. However, over 20 subsequent observations conducted both at 
the Parkes and the Lovell telescope between MID 56408 and 56675 
spending over 10 hr at the same position failed to re-detect the pul¬ 
sar. PSR 11759—24 was detected again since MID 56728 with both 
the Lovell and the Parkes telescopes and remains visible at the time 
of writing. 


Given the high DM of 772(14) cm"® pc of PSR 11759-24, 
it is highly unlikely that scintillation could be responsible for the 
intermittent detections. A tentative timing solution indicates that 
PSR 11759—24 could he in an eclipsing binary system with a long 
orbital period of the order of years (see Fig. [Jo). Similar exam¬ 
ples include PSR 11638—4725 which has a Porb of 5.3 yr and a 
~ 20 Mq ms companion, where the pulsar is undetectable in the 
radio wavelength for ~ 1 yr around periastron (Lyne et al., in prep); 
or PSR B1259—63 which has a Porb of 3.4yr and a ~ 10 Mq MS 
companion, where t he pulsar is undetect able during a 40-d eclipse 
behind the MS star dlohnston et al.|[l99^ . Further observations are 
needed to determine the nature of this pulsar and the cause for the 
intermittency. 


6.6 PSR J1847—0427, a pulsar with an extremely wide pulse 

PSR 11847—0427 has a very wide profile with emission over al¬ 
most the entire 360° of pulse phase. Given that most radio pul¬ 
sars have narrow emission patterns, PSR 11847—0427 is likely 
an aligned rotator with an angle between the magnetic and rota- 
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Table 6. Discovery parameters of 28 newly-discovered pulsars from the HTRU Galactic plane survey, which have not yet enough TOAs to produce a coherent 
timing solution. All of these pulsars have been assigned a temporary name containing only two digits of declination. We list the Galactic latitude (Z) and 
longitude (b), the sp in period (P) and the D M of these pulsars. Values in parentheses are the nominal Icr uncertainties in the last digits. The distances are 
derived according to lCordes & Lazid j2002h . 


PSR name 

R.A. (J2000) 

(hm) 

Dec. (J2000) 
(°') 

1 

(°) 

b 

(°) 

P 

(ms) 

DM 

(cm“^pc) 

Dist 

(kpc) 

J1227-63 

12:27.2(5) 

-63:09(7) 

300.20 

-0.41 

444.57796(6) 

359(4) 

8.4 

J1255-62 

12:55.3(5) 

-62:48(7) 

303.37 

0.06 

182.91150816(3) 

713(4) 

13.2 

J1349-63 

13:49.3(5) 

-63:56(7) 

309.24 

-1.77 

373.0340(1) 

478(3) 

9.6 

J1532-56 

15:32.3(5) 

-56:32(7) 

324.04 

0.01 

522.97708597(13) 

282 

4.3 

J1612-49 

16:12.9(5) 

-49:27(7) 

332.99 

1.28 

192.68718417(14) 

439(5) 

7.9 

J1612-55 

16:12.1(5) 

-55:09(7) 

328.99 

-2.78 

846.907(3) 

312(8) 

5.7 

J1627-49 

16:27.7(5) 

-49:54(7) 

334.37 

-0.71 

623.678(3) 

594(16) 

6.8 

J1627-51 

16:27.5(5) 

-51:08(7) 

333.46 

-1.54 

439.6840(5) 

201(4) 

3.5 

j 1634-49 

16:34.5(5) 

-49:52(7) 

335.13 

-1.49 

684.93649236(2) 

649(12) 

8.8 

J1638-44 

16:38.3(5) 

-44:40(7) 

339.43 

1.54 

568.0566269(3) 

494(5) 

6.7 

J1658-47 

16:58.4(5) 

-47:12(7) 

339.74 

-2.77 

369.3649(14) 

533(5) 

11.3 

J1710-37 

17:10.5(5) 

-37:30(7) 

348.82 

1.33 

792.65922183(3) 

346(15) 

4.7 

J1718-41 

17:18.2(5) 

-41:07(7) 

346.74 

-1.98 

548.4192(12) 

354(7) 

5.4 

J1720-36 

17:20.9(5) 

-36:53(7) 

350.51 

0.02 

92.13212(3) 

341(1) 

4.2 

J1723-38 

17:23.0(5) 

-38:20(7) 

349.56 

-1.15 

765.79581473(10) 

472(2) 

5.8 

J1730-34 

17:30.1(5) 

-34:48(7) 

353.28 

-0.36 

99.82956(15) 

628(5) 

7.0 

J1732-35 

17:32.5(4) 

-35:05(7) 

353.31 

-0.94 

126.689998553(5) 

340(2) 

5.1 

J1741-34 

17:41.9(5) 

-34:19(7) 

355.00 

-2.15 

875.137(2) 

241(8) 

4.1 

J1743-35 

17:43.1(5) 

-35:32(7) 

354.09 

-3.00 

569.980(9) 

174(5) 

3.3 

J1746-27 

17:46.0(5) 

-27:51(7) 

0.97 

0.49 

487.52760873(8) 

422(9) 

5.2 

J1748-30 

17:48.0(5) 

-30:17(7) 

359.12 

-1.14 

382.7347(6) 

584(5) 

8.1 

J1756-25 

17:56.7(5) 

-25:28(7) 

4.24 

-0.34 

855.986(2) 

706(8) 

9.6 

J1759-24 

17:59.4(5) 

-24:02(7) 

5.79 

-0.16 

1457.739(11) 

772(14) 

10.9 

J1819-17 

18:19.5(5) 

-17:05(7) 

14.15 

-0.90 

2352.135(15) 

67(23) 

1.9 

J1830-10 

18:30.1(5) 

-10:39(7) 

21.04 

-0.17 

245.260(3) 

203(2) 

3.6 

J1834-09 

18:34.6(5) 

-09:15(7) 

22.79 

-0.51 

512.3273(10) 

529(7) 

6.4 

J1836-11 

18:36.1(5) 

-11:17(7) 

21.15 

-1.77 

496.1931(11) 

326(8) 

5.2 



Figure 10. A tentative orbital solution for PSR J1759—24, with spin period 
plotted against the date of the observation in MJD. 


tional axis, a ~ 0°. This would imply that we are sampling the 
emission pattern along a ring around the magnetic axis that is en¬ 
tirely within the beam of emission. On top of this, the profile of 
PSR J1847—0427 exhibits a second component located at a rotation 
phase ~ 180° from the main pulse (see r elevant panel in F ig. [^, 
similar to that observed in PSR B0826—34 ( iBiggs et aljl988h . This 
second component could be due to oblateness in the emission beam. 
If future observations show evidences of drifting subpulses like the 
case of PSR B0826—34, this could provide insight into the con¬ 
ical beam of emission. Further measurements of the polarisation 
properties and observations of any frequency evolution can help to 
solve for the emission geometry, and would also allow us to con¬ 


strain the angle a. These can potentially lead to studies of pulsar 
magnetosphere and testing of emission mechanism models. 


6.7 Two intermittent pulsars 

PSRs J1227-63 (P = 444.6 ms) and 11349-63 (P = 373.0 ms) 
were both discovered with constant radio emission over the en¬ 
tire course of the 72-min survey observations. Including the dis¬ 
covery, the subsequent confirmation and timing observations as 
well as a weak detection in PMPS archival data for the case of 
PSR J1349—63, these two pulsars have been observed a total of 6.9 
and 13.4 hr respectively. However, they are not detectable in some 
of the follow-up observations. 

Given the high DM of 359(4) and 478(3) cm~® pc of these two 
pulsars, it is highly unlikely that scintillation could be responsible 
for the non-detections. Instead, they are potentially nulling pulsars 
with long nulls lasting up to the ~l-hr integration length of the 
follow-up o bservations, in which the pulsar is completely switched 
off (see e.g.. lKjamer et al.ll2006al) . Using Equation (O and the ob¬ 
served pulsar duty cycle at discovery, we derive upper limits on the 
flux densities corresponding to the non-detections in the longest 
follow-up observations of PSRs 11227—63 and 11349—63 to be 
0.08 and 0.04 mjy, respectively. 

It is pre-mature to calculate the periodicities in the switch¬ 
ing behaviour, but based on current data, we estimate that 
PSRs J1227—63 and J1349—63 are respectively visible, on aver¬ 
age, for less than 39 and 20 per cent of the time at 1.4 GHz. Fur- 
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Distance (kpc) 

Figure 11. Luminosity versus distance of the 54 newly-discovered pulsars 
with flux density measurements (red stars) and the published pulsars taken 
from PSRCAT (black dots). 

ther observations are necessary to determine the nulling fraction of 
these intermittent pulsars. 

7 COMPARING WITH KNOWN PULSAR POPULATION 

7.1 Luminosity 

54 out of the 60 newly-discovered pulsars have calibrated flux den¬ 
sity measurements and we have inferred their luminosities as re¬ 
ported in column 8 of Table We caution that the derived lumi¬ 
nosities are dependent on our knowledge of the pulsar distances. In 
this work we have calcu lated pulsar distances b ased on the NE2001 
electron density model dCordes & Lazidl2002h . which is typically 
thought to have an associated uncertainty of 25 per cent for each 
DM-derived distance. Assuming that this uncertainty is of compa¬ 
rable magnitude irrespective of the line of sight, we compare our 
new luminosity derivations to the luminosity distribution of the 
known pulsars as listed on PSRCAT. It is interesting to see from 
Fig.HUthat the HTRU Galactic plane survey is indeed probing the 
low luminosity region particularly for distant pulsars, as expected 
from the improved frequency and time resolution as well as the 
longer integration time of 72 min. One exception is PSR11759—24 
which has the highest luminosity among the newly discovered pul¬ 
sar. The intermittent detections associated with the eclipsing region 
have likely prevented this system (and possibly other intermittent 
pulsars along the Galactic plane) from being discovered previously. 

7.2 Distance 

No nearby pulsars have been discovered thus far, out to a distance 
of roughly 2kpc (see Tables |4] and |6l(. This null result might seem 
surprising at first; plus a quick examination of Fig. [TT] would sug¬ 
gest new discoveries to populate the region between the blue and 
green lines of constant flux densities. However, consider that the 
sampled sky volume increases with D^, for a distance D. If we as¬ 
sume that pulsars are uniformly distributed in the Galaxy, there are 


thus more pulsars at larger distances. Furthermore, a uniform pulsar 
distribution is certainly not a good assumption. In fact towards the 
inner Galactic plane, the pulsar density should increase as our line- 
of-sight crosses more spiral arms and approaches the Galactic cen¬ 
tre. This uneven Galactic distribution of pulsars further contributes 
to a higher number pulsars at larger distances for the HTRU Galac¬ 
tic plane survey region. Indeed, for the 713 known pulsars with a 
published luminosity within the HTRU Galactic plane survey re¬ 
gion, only 27 (i.e., 3.8 per cent) lie within a distance of 2 kpc. A 
Kolmogorov-Smirnov (KS) test comparing the newly-discovered 
pulsar distances with that of known pulsars within the same sky 
region of the HTRU Galactic plane survey indicates that the two 
samples are consistent with being drawn from the same population, 
with a p-value of 0.4 (i.e., probability > 0.05). 

However, if the complete HTRU Galactic plane survey pro¬ 
duces a significantly smaller percentage of nearby discoveries than 
3.8 per cent, it could indicate that we have completed the nearby 
pulsar population, or at least reached a point where the yield of 
pulsar surveys are reducing as we are no longer flux limited. This 
has important implications to future pulsar surveys such as those 
to be conducted with MeerKAT and the SKA. Any pulsar surveys 
targeting the Galactic plane will have to go to higher observing fre¬ 
quency which could help with the discovery of pulsars that would 
be otherwise undetected due to scattering. A yet longer dwell time 
or a telescope with larger collecting area would thus be needed to 
compensate for the typically negative spectral index of pulsars as 
they get weaker at a higher observing frequency. 

7.3 Characteristic age 

Young pulsars (conventionally defined as having characteristic ages 
less than 100 kyr) are expected to be located not far from their birth 
places hence mostly populating the Galactic plane region. Indeed 
the previous Galactic plane survey, the PMPS, discovered a sam¬ 
ple of pulsars with an avera ge characteristic ag e lower than the 
previously-known population IMorris et alj2002) . and these PMPS 
discover ies account for abou t half of the young pulsars currently 
known l lKramer et alj l2003h . No young pulsa rs have been dis - 
covered in the HTRU medium-latitude survey dBates et al.ll2012h . 
This null result was attributed to the fact that the HTRU medium- 
latitude survey has a higher minimum detectable flux density of 
0.2 mi y compared to the PMPS of 0.15mJy along the Galactic 
plane, bates et ^ (1201 2 h predicted that young pulsar discoveries 
are instead expected from the HTRU Galactic plane survey, given 
the higher sensitivity of this project part. 

Of the 60 newly-discovered pulsars presented here, 32 now 
have coherent timing solutions obtained over a time span of one 
year or more, allowing us to infer the preliminary characteristic 
ages and to locate them on a P-P diagram as in Fig. [I2\ Note that 
we have not corrected for any potential contribution in the observed 
P from the Galac tic potential and f rom the transverse proper mo¬ 
tion of the pulsar dShklovskiill 19701) . T hese effects are expe cted to 
be significant only for recycled pulsars dCamilo et al.lll993) . 

There is a noticeable lack of young pulsars amongst the 32 
which have timing solutions as they all have Tc at least of the order 
of Myr (also see Table|4j. We caution that the inferred characteris¬ 
tic ages for recycled pulsars (shade d region in F ig. [12) a re known t o 
be unreliable estimates in general dTaurisl2012l ; lTauris et al.ll2012h . 
Hence, we disregard all recycled pulsars, both previously-known 
and newly-discovered from this survey (namely PSRs 11101—6424 
and J1244—6309), and we compute a KS test comparing the ages of 
the newly-discovered pulsars with that of the known pulsars within 
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Figure 12. P-P diagram showing the known pulsar population (black) 
within the HTRU Galactic plane survey region. The 32 new discoveries 
that have now preliminary timing solutions are denoted by red stars, except 
PSR J1829—1009 for which a 2a upper limit is shown by a blue star. Lines 
of constant characteristic ages are shown as red dotted lines, and it can be 
seen that the newly discovered pulsars appear to be of an older population. 
In addition, lines of constant spin-down energies are shwon as cyan solid 
lines and surface magnetic field strength as blue dashed lines. The death 
line as described in IChen & RudermanI il993ll is shown as green dot-dashed 
line. We define ‘recycled pulsars’ to be those with Sgurf ^ 3 X 10^® G as 
indicated by the shaded region. 

the same sky region of the HTRU Galactic plane survey. The re¬ 
sult indicates that the two samples are not consistent with being 
drawn from the same population, with a p-value of 0.005 (i.e., 
probability < 0.05). This is at odds with the expectation mentioned 
above. The identity of this older and less luminous population is in¬ 
triguing, and may account for the missing population of disrupted 
mildly-recycled puls ars along the Galactic plan e as predicted by 
llxmmer et all ^2004) and further investigated by iBelczvnski et al.l 
1 201 oh . a population that has as yet largely eluded pulsar surveys 
due to their small average flux density. We note that the Galactic 
plane location of this old population also implies that these pulsars 
have likely received small kicks at birth. 

In addition, there is a long standing debate regarding the re¬ 
lationship of radio luminosity of pulsars and their ages. In princi¬ 
ple, the radio emission only contribut es to a small portion of the 
energy budget of a pulsar. However, iNaravan & Ostrike3 l ll990h 
suggested that radio luminosity is proportional to the cube root 
of the spin-down luminosity of a pulsar (i.e., L ex P^^^/P). In 
other words, older p ulsars would have a lower radio luminosity. 
lUorimer et al, Il ll993l) sh owed that there is a l arge s catter in this re¬ 
lationship. Nonetheless. [Arzoumanian et alj ( l2002h argued that the 
spread is simply caused by viewing geometry and intrinsically there 
is a strong luminosity-age relationship that scales with the voltage 
drop above the polar cap. If this is true, it would provide an ex¬ 
planation for the discovery of these older pulsars, as the HTRU 


Galactic plane survey probes deeper in the low luminosity end of 
the Galactic pulsar population. 

Timing solutions for the remaining 28 pulsars currently with¬ 
out a P measurement as well as any future discoveries from the 
HTRU Galactic plane survey will be crucial to differentiate if the 
current lack of young pulsars is purely due to small number statis¬ 
tics. The study of the age distribution of this less-luminous pul¬ 
sar population probed by the HTRU Galactic plane survey should 
help to resolve some of the degeneracy of the arguments introduced 
above. 

8 A COMPARISON WITH THE ESTIMATED 

DETECTIONS 

Usi ng the PSRPOlP^ software based on t he pul sar population model 
by lUorimer et al.l J2()0^ . iKeith et ^ ( 1201 Oh estimated that the 
HTRU Galactic plane survey should detect 957 normal pulsar (here 
defined to be of spin period longer than 30 ms). We have conducted 
a similar simulation u sing the mo r e up-t o-date PsrPopPvF^ soft¬ 
ware as presented in [Bates et ^ ([20141) . By assuming the s ame 
power-law luminosity model as employed in [Keith et al.[([2010[) . we 
find that the expected number of norma l pulsar detections is 1020, 
statistically consistent with the result in [Keith et alj ( [201(lh . 

Scaling these estimated numbers down to reflect the current 
HTRU Galactic plane survey data processing progress of 50 per 
cent, it can be seen that of the order of 470 to 510 normal pul¬ 
sar detections are predicted. We note that, strictly speaking, such 
intrapolation is inaccurate, as it assumes pulsars to be uniformly 
distributed in the survey region. This is not true; for example, a 
higher density of pulsars should be observed along the line of sights 
that cut through any spiral arms. However, as shown in Fig.[T] this 
50 per cent of processed observations are drawn from random se¬ 
lection, hence this intrapolated estimation should hold. 

Among the 435 known pulsar re-detections, 427 are normal 
pulsars. Adding this to the 58 normal pulsars discovered thus far, 
we have a current normal pulsar detection rate of 485. The close 
match between the estimated and observed yield is satisfying, and 
would imply that our survey is performing as expected. In addi¬ 
tion, the ability of these simulation tools to closely reproduce the 
observed yield could be used as an indication of a sound under¬ 
standing of the underlying pulsar population parameters. Nonethe¬ 
less, we note that a more mature judgment can be formed when 
the complete HTRU Galactic plane survey is analysed in the near 
future. 

The detectability of the fastest-spinning MSPs (here defined to 
be of spin period ^ 30 ms) is subjected to different pulsar-searching 
selection bias and is often considered separately. PSRPOP and PSR- 
PopPy have predicted 51 and 43 MSP detections respectively for 
the complete HTR U Galactic plane su rvey. A dedicated study of 
MSP luminosity by [Levin et ^ ( [2013h suggested an even higher 
MSP yield of 68. Hence within the current 50 per cent of processed 
observations, these simulations predict an MSP detection rate be¬ 
tween 25 and 34. Among the 435 re-detected known pulsars, eight 
are attributed to previously known MSPs, and two additional MSP 
detections come from the new discoveries of PSRs J1101—6424 
(Section [6.lb and 11757—27 (Section [6.4b . These sum to a current 
total MSP yield of ten. It thus appears that we are roughly a factor 
of three short of MSP detections. 

http://psrpop.sourceforge.net 

https://github.com/samb8s/PsrPopPy 
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Perhaps the most relevant explanation here, is that we have 
only processed the HTRU Galactic plane data in a partially- 
coherent manner. The segmentation scheme means that we have not 
make use of the full sensitivity achievable by the complete 4300 s 
observations. In fact, running the MSP simulation again with PSR- 
PopPy for each of the parallel searches of halved, quartered, and 
one-eighth segment length, the corresponding predicted MSP de¬ 
tections would be reduced to 27, 22, and 16 respectively, for the 
entire Galactic plane sub-survey. The current number of ten MSP 
detections in 50 per cent of processed data then become at least 
consistent. As discussed in Section [3.2.3l coherently analysing the 
full 4300 s observation with acceleration search is our priority in 
future data reprocessing. 

Apart from that, a major caveat in the estimation of MSP de¬ 
tection is that, our understanding of the underlying MSP distribu¬ 
tion is still inadequate and hampered by the small number statistics 
of known Galactic MSPs. For this reason, MSP po pulation sim- 
ulation was no t considered in iLorimer et al.l ( l2006h . In addition, 
iLorimeil ( l20I3h pointed out potential issues involv ed in the MSP 
period model devised bv ICordes & Chem ofll 1 19 9^ whic h is used 
as an input parameter in PsrPopPy. iLevin et alj(l2013h found a 
steeper luminosity distribution and suggested a larger ^-he i ght fo r 
MSPs than what was used in the simulation by I Keith et al.l d201(t) . 
These disagreements in the input parameters mean that we should 
treat any predictions from MSP population simulations with care. 

Despite these uncertainties in the MSP population model, sev- 
eral other reasons co uld come into play reducing our MSP yield. 
lEatough et al.l ( 1201 3h have provided a detailed discussion on the 
possible causes of MSP non-detection. For instance, they suggested 
that the influence of unfavourable orbital phases where the acceler¬ 
ation was not constant, the emission beam having precessed out of 
our line of sight, intermittency, as well as human inspection error 
during pulsar candidate selection could have accounted for some of 
the missed MSPs. 

Furthermore, scattering broadening caused by the interstellar 
medium is particularly harmful for MSP detection. Any broadening 
of a time scale of ms already represents a significant portion of the 
spin period of an MSP, and the fact that MSPs tend to have larger 
duty cycles (see Table [TJ makes them more susceptible to profile 
broadening. Scattering is not well quantified and if this effect is 
underestimated, could partially explain the apparent low yield of 
MSPs. We expect that the large sample of pulsars (both normal 
and MSPs) to be detected by the HTRU Galactic plane survey will 
provide further information for a better quantification of the extent 
of scattering. 

The presence of RFI could also have reduced our sensitivity 
which would have particularly affected the detection of low lumi¬ 
nosity pulsars. However, as demonstrated by the close match of 
yield in normal pulsar detections, the HTRU Galactic plane survey 
does not seem to have any obvious general loss of sensitivity, hence 
RFI should not be a significant factor for the explanation of our lack 
of MSP detections. 


9 DISCUSSION AND CONCLUSION 

The HTRU survey is a blind puls ar survey of the so uthern sky with 
the 64-m Parkes radio telescope dKeith et al.ll2010l) and a twin sur- 
vey of the northe rn sky with the 100-m Effelsberg radio telescope 
dBarr et alj|2013l) . The HTRU surveys have benefited from recent 
advancements in technology and provide unprecedented time and 
frequency resolution, making the HTRU the first complete all-sky 


survey and the best pre-SKA survey. Thus far, over 170 newly- 
discovered pulsars have been reported from the HTRU survey 
with Parkes dKeith et atese^ DI 2 OI l|.lKeith et al.ll20li 


Bates et aDl2012l ; Burgav et^ 20131 : Ng et al. 12014 : Bates et aH 

2015h . induding 33 MSPs. 

The HTRU low-latitude Galactic plane survey promises to 
provide the deepest large-scale search conducted thus far for the 
Galactic plane regi on, where the most rel ativistic binaries are ex¬ 
pected to be found dBelczvnski et alj2002l) . In addition, this survey 
will represent a unique record of the Galactic plane with high yet 
uniform sensitivity, enabling an unbiased Galactic census to ex¬ 
plore the true boundaries of pulsar phase space. The discoveries 
from the HTRU will also provide valuable knowledge of the Galac¬ 
tic pulsar population for the planning of survey strategies with the 
SKA. 


Searching for pulsars in high resolution observations with long 
integration lengths is computationally intensive, and the depth to 
which the data can be explored is limited by the available com¬ 
puting resources. In order to improve the searching algorithm, we 
introduce two REI mitigation techniques in the time and Fourier 
domain. The Fourier domain method has been shown to be effec¬ 
tive and for some cases pulsars have been discovered only when 
this technique is applied. Furthermore, we present the implementa¬ 
tion of a novel partially-coherent segmented acceleration search al¬ 
gorithm, based on the time-domain re-sampling method. This seg¬ 
mented search technique aims to increase our chances of discov¬ 
eries of highly-accelerated relativistic short-orbit binary systems. 
We show that a rorb ~ 0.1 results in the highest effectiveness of 
the constant acceleration approximation, hence we split the long 
data sets of the HTRU Galactic plane survey into segments of var¬ 
ious lengths and search for the binary systems using an accelera¬ 
tion range appropriate for the length of the segment. We push the 
maximum achievable acceleration value to 1200 ms“^. Within this 
range we can find a double-pulsar-like system deeper in the Galaxy 
with orbital periods as short as 1.5 hr, and explore the parameter 
space occupied by expected pulsar-black hole systems. Of the or¬ 
der of 1.8 X 10^"^ computational operations per data set are required 
by the search configuration adopted in this survey. The high com¬ 
putational requirements of such data-intensive astronomy also act 
as a test bed for SKA technologies. 

Analysis of 50 per cent of the HTRU Galactic plane obser¬ 
vations with the partially-coherent segmented acceleration search 
outlined above has resulted in 912 re-detections of 435 previously- 
known pulsars, demonstrating that the survey is performing as ex¬ 
pected. In addition, we present the discovery of 60 pulsars, of which 
two are fast spinning pulsars with periods less than 30 ms. One of 
the two pulsars, PSR J1101—6424, is likely a descendant from an 
IMXB. Its fast spin period of 5.1 ms indicates a rare example of 
full-recycling from a long mass-transfer phase, making it only the 
second known IMXB which has evolved from a Case A Roche lobe 
overflow. The second one, PSR J1757—27, is likely to be an iso¬ 
lated pulsar. The most adopted formation scenario of an isolated 
MSP requires the pulsar to have a very fast spin period of the or¬ 
der of a few ms, at odds with the unexpectedly long spin period of 
17.7 ms in the case of PSR 11757—27. Furthermore, we report on 
the discovery of PSR 11244—6359, a mildly-recycled 147.3 ms pul¬ 
sar in a 17.2-d binary orbit, bridging the normal and the recycled 
pulsar population on a P-P diagram. We report on the eccentric 
(e = 0.09) binary pulsar PSR J1755—25 with a heavy companion. 
We also report on the discovery of PSR J1759—24. This pulsar is 
likely to be in an eclipsing binary system with a long orbital pe¬ 
riod of the order of tens of years. We report on the discovery of 
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an aligned rotator PSR J1847—0427 whose pulse profile contains 
a main pulse and a second component separated by 180°, together 
forming a wide profile with emission over almost the entire 360° of 
longitude. We also report on the discovery of two intermittent pul¬ 
sars, PSRs J1227—63 and J1349—63. From the current discovery 
rate, extrapolation shows that the Galactic plane survey will result 
in roughly a further 60 discoveries. 

We note that the sky region of the HTRU Galactic plane sur¬ 
vey has a complete overlap with that of the PMPS, as well as the 
HTRU medium-latitude survey. Despite the common sky coverage 
and the high success rates of these previous pulsar surveys, pulsar 
discoveries are still continuously being made from this survey. The 
improved dynamic range of the BPSR digital backend over the pre¬ 
vious analogue filterbank system employed by the PMPS, as well 
as the long integration length of 72 min of this survey, both account 
for the high survey sensitivity crucial for the discoveries of many 
less luminous pulsars. Nonetheless, some of the newly-discovered 
pulsars presented in this paper are later found to be detectable in 
the archival data. The discoveries of these relatively bright pulsars 
which have eluded previous survey efforts could be attributed to 
advancements in the data processing techniques. Particularly, the 
improvement in RFI mitigation techniques have reduced the num¬ 
ber of false-positive candidates. Moreover, the ‘multiple-pass’ na¬ 
ture of the partially-coherent segmented acceleration search has in¬ 
creased the chance of discovering pulsars by avoiding parts of the 
observation in which they might be less detectable. In principle, 
the acceleration search would also provide sensitivity to highly- 
accelerated relativistic pulsar binaries. However, none of the newly- 
discovered pulsars reported here was found at an orbital phase with 
noticeable orbital acceleration. 

Of the newly-discovered pulsars presented here, 32 now have 
coherent timing solutions, which allows us to infer their prelimi¬ 
nary characteristic ages. There is a noticeable lack of young pul¬ 
sars among these pulsars as they all have Tc at least of the order 
of Myr. Using a KS test to compare the characteristic ages of this 
newly-discovered population with that previously known, we show 
that they are inconsistent of being drawn from the same population. 
The old age and low E/d? of these pulsars make Fermi associations 
unlikely. Timing solutions for these pulsars and any future discov¬ 
eries from the HTRU Galactic plane survey will be key to study the 
characteristic age distribution of this less-luminous pulsar popula¬ 
tion and might help to resolve some of the long standing arguments 
regarding the relationship of radio luminosity of pulsars and their 
ages. 

A comparison with the estimated survey yield shows that we 
currently have a close match between the estimated and observed 
yield of ‘normal’ pulsars. This is satisfying and would imply that 
our survey is performing as expected. We appear to be roughly a 
factor three short of MSP detections, and in Section [8] we explore 
the potential causes of the missing MSPs. We point out that the 
most relevant explanation is the fact that we have only processed 
the HTRU Galactic plane data in a partially coherent manner. The 
segmentation scheme means that we have not make use of the full 
sensitivity of the survey. Coherently analysing the entire observa¬ 
tion length with an appropriate acceleration range is thus our first 
priority in future data reprocessing. 
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APPENDIX A: TECHNICAL DETAILS OE THE 
SEGMENTED ACCELERATION SEARCH PIPELINE 

AI Determining the optimal Torb 

To quantify the effectiveness of a constant acceleration approxi¬ 
mation versus varying roib, we employed two observations of the 
double pulsar system PSR J0737—3039A as test data sets. The dou¬ 
ble pulsar is the most relativistic pulsar binary system known, with 
very high maximum orbital acceleration of the order of 250 ms“^ 
and a short Porb of 2.45 hr. The orbital eccentricity is reasonably 
small (e = 0.088), so it closely reproduces the simplest orbital 
motion; a sinusoidal v{t) of a circular orbit. Two test observations 
of the double pulsar system were carried out at Parkes. One was 
collected with the same observational set-up as the HTRU Galactic 
plane survey, where fsamp = 64 ps (hereafter test d ata set 1), and 
the oth er is identical to the test data set presented in lEatough et al.l 
(l2013h . where fsamp = 80 ps (hereafter test data set 2). The two 
different fsamp allow different rorb to be probed, as we segment 
the test data sets into increasingly shorter tint, each with tifft of 
2*^ to allow for maximum computational efficiency of the FFTs, 
where k = 25, 24, 23, 22, 21, 20. 

For each rorb, we conduct a series of acceleration searches on 
a subset of the observation with the relevant length, incrementing 
at every 50 s across the observed orbital phase. In Figs. lAll and lA2l 
we plot the detected spectral S/hQ of the double pulsar in black 
and the recovered orbital acceleration in other colours. Panels (a) 
of both figures show the shortest rorb for each respective test data 
set. From the red data points it can be seen that these short seg¬ 
ments with small rorb each contain only a tiny fraction of the or¬ 
bital motion of the double pulsar. Hence the analyses are not very 
sensitive to the trial acceleration value, leading to large fluctuations 
in the recovered orbital acceleration which does not follow the pre¬ 
dicted grey curve closely. Nonetheless, the double pulsar has been 
detected throughout the orbital phase with roughly consistent spec¬ 
tral S/N, an indication that the constant acceleration approxima¬ 
tion has been equally effective irrespective of orbital phase. As the 
coherent segment length gets progressively longer, as shown from 
panels (b) to (d), the acceleration searches become more success¬ 
ful in recovering the predicted orbital acceleration and the detected 
S/N also improves by roughly v/2 as expected by the radiometer 
equation. Panels (e) represent yet longer values of rorb which ex¬ 
ceed 0.1. At phases 0, 1, and 0.5, where the orbital acceleration is 
closest to being constant (i.e., the acceleration derivative d « 0), 
improvements in the detected S/N are still observed compared to 
the previous rorb ■ However, the increase is less thanv/2, reflecting 


The .spectral S/N plotted here is the best S/N from the summing of the 
five harmonics mentioned in Section[3] 
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Figure Al. Detected orbital acceleration and S/N at various orbital phases 
for test data set 1. Each panel corresponds to progressively longer values of 
r-„b, from (a) 0.0076, (b) 0.015, (c) 0.030, (d) 0.061, (e) 0.12 to (f) 0.24. 
The black bar on the bottom left of each panel depicts the length of each 
segment to be searched coherently. We slide this search window across the 
test data set incremented by every 50 s, and we plot the acceleration search 
result at the middle of each segment. The recovered orbital acceleration 
values are plotted with colours, and the detected S/N is plotted in black. 
The smooth grey line shows the predicted acceleration curve. 

that the rorb is now hecoming too large resulting in spectral smear¬ 
ing which in turn reduces the detected spectral S/N, making the 
constant acceleration approximation less effective. At phases 0.25 
and 0.75, where the orbital acceleration is increasing the most (i.e., 
a significant a), the detected S/N is worse than the shorter rorb - Fi¬ 
nally, panels (f) show the longest rorb where the degradation due to 
spectral smearing out-weighs the gain in S/N due to longer coher¬ 
ent segments, producing lower S/Ns at all orbital phases. The dras¬ 
tic drop in S/N immediately away from the orbital phases where 
d ~ 0 is particularly noticeab le. These are in agreement with the 
results of lEatough et alj ( l2013l) . 

A2 Computational power considerations 

For this partially-coherent acceleration search, we adopt the same 
accele r ation step size as descr ibed by Equation (3) of lEatough et al.l 
( l2013h . lEatough et alj (12013h pointed out that pulsars with longer 
spin periods are less susceptible to a wrong acceleration trial (see 
Fig. 3 in their publication), and hence this acceleration step size 
is effectively slightly oversampling for these slow spinning bina¬ 
ries. In turn, the number of computational operations required for 
Fourier transforming each real time series, Cs, can be quantified as 

C7s OC An X fsamp X Ttsamp X log 2 (ttsamp ) • (Al) 


Figure A2. Detected orbital acceleration and S/N at various orbital phases 
for test data set 2. Each panel corresponds to progressively longer values of 
r„b, from (a) 0.0095, (b) 0.019, (c) 0.038, (d) 0.076, (e) 0.15 to (f) 0.30. 

Note that because we progressively downsample in time after 
2x the diagonal DM is reached (see e.g.. Section 6.1.1.2 of 
iLorimer & Kramen l2005h . fsamp does not take the same value 
across different DM values, as indicated in Fig. Due to con¬ 
straints of computational resources, we further downsample all ob¬ 
servations to a fsamp ^ 256 qs. This might reduce our detectabil¬ 
ity towards MSPs with spin periods ^ 1 ms, but st i ll com pares 
favourably to the fsamp = 1 ms used in lEatough et al.l ( l2013ll . 

The total computational operations required to Fourier trans¬ 
form for each data set, Ctot, then becomes 

s 

Ctot OC ^ Cs X trialsDM x s, (A2) 

s=l,2,4,8 

where trialsoM is the number of DM trials and its value is indicated 
in Fig. [6] Summing the four configurations with segments of s = 
1, 2 , 4,8, we have of the order of Ctot ~ 1-8 x 10^^ operations. 
Such processing typically requires ~620 CPU core hours on a sin¬ 
gle Intel Xeon Sandy Bridge node computer, for the analysis of one 
beam of an observation. Fig. [A^ illustrates the relative processing 
time needed for each of the four configurations as a function of the 
number of DM trials. 


APPENDIX B: CANDIDATE CONEIRMATION AND 
GRIDDING STRATEGY 

When a promising pulsar candidate is identified, it is necessary to 
conduct a confirmation observation at the telescope to verify if the 
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Number of DM trials 


Figure A3. The relative processing time required for each of the four con¬ 
figurations (s = 1, 2,4, 8) parallel searches, normalised by the total pro¬ 
cessing time of all segments together. We plot these curves as a function of 
the number of increasing DM trials, and we show the corresponding DM 
values. The relevant fsamp is also marked. 


candidate is a genuine pulsar. A successful re-detection confirms 
that the pulsar can he seen within the Gaussian beam of the re¬ 
ceiver which has a FWHM of 14'.4 at Parkes. However, a better 
localisation of the pulsar position is desirable. Pulsars north of dec¬ 
lination —35° are followed-up with a timing campaign at the Lovell 
TelesconJ^ 

Major pulsar surveys at the Parkes telescope, such as the 
PMPS as well as the medium- and high-latitude parts of the 
HTRU, carried out a ‘N-E-S-W’ gridding strategy as described in 
iMorris et alj ( l2002h for achieving a better positioning of any newly- 
discovered pulsar. In addition to one re-observation at the exact po¬ 
sition of the discovery (D), a grid surrounding the discovery posi¬ 
tion is carried out. An offset of = 9' is applied in each direc¬ 
tion of North (N), East (E), South (S), West (W) from the discov¬ 
ery position. This offset has been chosen so that these four grids 
form a tight square pattern through the discovery position fFie. lBlI 
upper panel). The Gaussian beam of the receiver would overlap 
enough between each grid position to ensure that the uncertainty 
of the pulsar position is less than a single beam width. The respec¬ 
tive detected S/N from the ‘N-E-S-W’ pointings can then be used 
to estimate the true position of the pulsar. 

The long integration length of the HTRU Galactic plane ob¬ 
servations implies that, adopting this ‘N-E-S-W’ gridding scheme 
as it is used in the other parts of the HTRU survey, a total of 
72 min x 5 = 6 hr would be needed to conhrm each pulsar candi¬ 
date which is highly inefficient. We hence employed an optimised 
strategy with a ‘Ring-of-3’ grid (Eig. lBlI lower panel) to minimise 
telescope time. Although this has a slightly worse sky coverage 
as compared to the ‘N-E-S-W’ scheme, the factor of 5/3 gain in 
time is desirable. Furthermore, the integration time for each con- 


In practice, for the declination between —30° and —35°, only bright 
pulsar.s are followed-up at the Lovell Telescope. This is because of the short 
visible hours, as well as the challenges associated with such low elevation 
observations, namely the need for low wind conditions and the potential 
high spillover, (see Section |3, which has a smaller FWHM of 12'. In ad¬ 
dition, an accurate position ensures efficient timing observations, as this 
maximises the S/N and hence reduces the telescope time required. 
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Figure Bl. Illustration of the gridding configuration. The light grey circles 
indicate beam pattern of the 13-beam Multibeam receiver with a FWHM 
of 14'.4. The upp er panel shows the gridding strategy as described in 
iMorris et alJ J2002h where five positions (‘N-E-S-W’ and the discovery po¬ 
sition ‘D’) are required each offsetting the discovery position by (5i = 9'. 
The lower panel shows the ‘Ring-of-3’ configuration used in the gridding in 
this survey, where only three positions (A-B-C) are required each offsetting 
the discovery position by S 2 = l/ZS X FWHM Ri 8'.3. 


firmation observation is scaled down from the nominal discovery 
S/N to achieve an expected S/N of 10 in the confirmation. Archival 
data (including overlapping observations from the HTRU medium- 
latitude survey and the PMPS observations) are checked before any 
gridding is carried out, in order to improve a priori knowledge of 
the true position of the pulsar. In practice, we find that almost all of 
the 60 newly-discovered pulsars were successfully confirmed fol¬ 
lowing this ‘Ring-of-3’ scheme, affirming the efficiency over pre¬ 
vious gridding strategies. 


APPENDIX C: PREVIOUSLY-KNOWN PULSARS WITHIN 
THE SURVEY REGION 

Cl A database of all re-detections 

In all, 912 re-detections of 435 previously known pulsars were ob¬ 
served in the HTRU Galactic plane survey data. Their details are 
listed in Table lClI ffull version available as Supporting Information 
with the online version of the paper). 
































HTRU—XII. The Galactic plane survey 23 


Table Cl. The 435 previously known pulsars that have been re-detected in the 50 per cent processed data of the HTRU Galactic plane survey. The file ID 
(which stands for the date and the time in UT of the observation) for which the pulsar was re-detected is listed, as well as the Galactic longitude (Z) and latitude 
(b) that corresponds to the central position of the receiver beam of that observation. We list also the offset between this and the true position of the pulsar. Note 
that in the case when a pulsar was re-detected multiple times in different observations, we have listed only the re-detection which was closest to the pulsar, i.e. 
that with the smallest offset. We list the observed spin period (P) and the DM, as well as the expected flux density at 1.4 GHz (Sexp) and the expected S/N 
(S/Nexp) calculated according to the radiometer equation. For the pulsars where no flux density is published on PSRCAT, we indicate with an asterisk (*). 
Finally, we list the observed S/N (S/Notis)- The folllowing is a sample of the full table, which is available as Supporting Information with the online version 
of the article. 


PSR name 

file 

Pointing/beam 

Z 

n 

b 

(°) 

offset 

n 

-^obs 

(ms) 

(cm~^pc) 

‘S'exp 

(mJy) 

S/Nexp 

S/Noba 

B0959-54 

2011-10-11-18:36:49/01 

280.230 

0.060 

0.03 

1436.619* 

129.3 

6.10 

2100.0 

1384.1 

BlOll-58 

2011-01-29-15:34:44/08 

283.587 

-2.193 

0.13 

819.922 

374.8 

0.73 

184.0 

129.8 

B1014-53 

2010-12-08-17:01:17/08 

281.170 

2.454 

0.03 

769.584 

65.5 

0.78 

204.0 

81.4 

B1015-56 

2011-04-25-06:10:43/13 

282.712 

0.237 

0.11 

503.462 

434.5 

1.64 

378.0 

217.1 

B1030-58 

2011-05-19-06:55:11/05 

285.963 

-0.940 

0.07 

464.210 

420.1 

0.73 

187.0 

154.6 

B1036-58 

2011-07-14-00:46:57/01 

286.300 

0.060 

0.08 

661.994 

70.3 

0.55 

123.0 

264.6 

B1039-55 

2011-07-07-02:30:49/03 

285.235 

3.100 

0.11 

1170.865 

299.4 

0.34 

75.5 

160.2 

B1044-57 

2013-02-25-12:20:14/12 

287.153 

0.683 

0.10 

369.431*1 

211.6 

0.65 

115.0 

44.1 

B1046-58 

2012-08-07-00:41:51/03 

287.335 

0.660 

0.12 

123.714 

128.8 

3.03 

462.0 

531.1 

B1054-62 

2011-07-18-00:07:15/12 

290.413 

-2.967 

0.12 

422.450 

322.2 

9.95 

1250.0 

1319.6 


1 Spin period originally detected at the third harmonics, i.e. at one third of the fundamental spin period listed here. 
Spin period originally detected at the half harmonics, i.e. at twice the fundamental spin period listed here. 


C2 Non-detections 

Apart from the magnetar PSR J1622—4950 discussed in Section|4l 
three further non-detections are related to a transient anomalous X- 
ray pulsar (AXP), probably in a radio quiet state. The non-detection 
of PSR J1747—2958 can be explained by the high flux variability 
due to scintillation by the associated nebula. PSR J1746—2849 is 
known to have a very broad profile and the fact that it is within 
1° from Sgr A* explains its significant scattering time scale of 
Tsc,i.4GHz > 226 ms, potentially preventing it from being de¬ 
tected in our survey. A further four non-detections are associated 
with pulsars with known intermittency and we might have observed 
them at a time when they were not emitting. 

Excluding the magnetar, the AXP, the pulsars in local vari¬ 
able environment and the intermittent pulsars, twelve further non¬ 
detections are listed in the last three panels of Table IC2| Five non¬ 
detections are attributed to the presence of RFI in individual sur¬ 
vey observations or confusion with red noise. For PSRs J1644—44 
and J1644—46, the published declinations are not well-determined 
(E. Keane; private communication). This means that the actual 
S/Nexp might be lower if the true offset between the pulsar and 
the observed position is larger. The pulsars listed in the last panel 
of Table IC2l have unexpectedly weak S/Ns. Their relatively high 
DMs imply that scintillation is unlikely to have an influence on 
the detectability. It is not clear whether this is a result of an over¬ 
estimation in their catalogue flux densities or if some unknown phe¬ 
nomena such as nulling has occurred intrinsic to the pulsar. For 
PSR J1746—2850, other observers have also failed in re-detecting 
any emission from other telescopes (G. Desvignes, P. Lazarus; pri¬ 
vate communication). 

C3 Re-detections of previously-known binary pulsars 

Table lC3l lists all previously-known binary pulsars re-detected thus 
far. Two known binary systems (namely PSRs J1525—5545 and 
J1807—2459A) which were undetected in our initial ‘standard 
search’ pipeline were later found from the ‘partially-coherent seg¬ 
mented acceleration search’ pipeline, and some other binaries have 
been detected with much higher S/Ns (for example in the case of 


PSR J0737—3039A) from the acceleration search. These examples 
illustrate that the acceleration search algorithm allows us to detect 
fast binary systems which we would be insensitive to otherwise. 
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Table C2. Previously known pulsars with an S/Nexp > 9 that have been missed in the 50 per cent processed observations of the HTRU Galactic plane survey. 
We list the file ID (which stands for the date and the time in UT) of the closest HTRU observation together with the positional offset from the catalogue pulsar 
position. The catalogue spin period and DM of the missed pulsars are also shown, together with the expected S/N (S/Nexp) as well as the recovered S/N when 
folding observation with catalogue ephemerides (S/Neph)- 


PSR name 

Observation ID 

offset 

P 

DM 

S/Nexp 

S/Ne 

ph Comments 


(Pointing/Beam) 

(“) 

(s) 

(cm“® 

pc) 


J1622-4950 

2010-12-29-03:31:28/08 

0.19 

4.326100 

820.0 

20.2 

— 

Radio magnetar (Levin et al.ll201C|) 

J1809-1943 

2011-05-17-17:28:45/10 

0.21 

5.540354 

178.0 

88.4 

- 

AXPXTF,.11810-197 jCamilo et nlj20n6|l 


2011-06-26-12:10:21/03 

0.22 



78.1 

- 



2011-07-02-12:25:09/02 

0.14 



330.1 

— 


J1747-2958 

2013-01-07-03:05:18/08 

0.12 

0.098814 

101.5 

12.8 

— 

Associated with the ‘Mouse’ radio nebula. Flux variability due to 
interstellar scintillation (Camilo et alj2002tl 

J1746-2849 

2011-06-26-10:07:53/08 

0.10 

1.478480 

1456.0 

20.5 

— 

Very broad profile with r^r ^ 4 oH’?, > 226 ms and proximity to 
Sgr A* iDeneva et al J2009tl 

J1634-5107 

2011-05-07-19:53:44/13 

0.07 

0.507356 

372.8 

46.4 

- 

Known inteiTnittencv (O’Brien et all2006l) 

J1726-31 

2011-05-08-17:04:49/03 

0.042 

0.123470 

264.4 

12.2 

— 

Known intermittency, detection probabilitv 20% (Knispel et alj 
|2013b, published declination not well-determined. 

J1808-1517 

2013-04-02-18:57:44/10 

0.13 

0.544549 

205.0 

20.8 

— 

Known intermittencv lEatouah et al.l2013l) 

Known nulling. (Wana et alJ2007h 

B1713-40 

2013-04-01-14:21:48/05 

0.23 

0.887710 

308.5 

715.0 

- 

B1727-33 

2013-02-01-00:30:51/13 

0.21 

0.139460 

259.0 

36.6 

- 

Observation badly affected by RFI. 

J1730-3353 

2013-02-01-00:30:51/13 

0.12 

3.270242 

256.0 

38.4 

- 

Observation badly affected by RFI. 

J1845-0316 

2011-04-25-19:36:08/03 

0.13 

0.207636 

500.0 

20.9 

13.2 

Observation affected by RFI. 

J1840-0840 

2013-01-06-02:28:11/12 

0.21 

5.309377 

272.0 

16.5 

24.1 

Long period pulsar, confusion with red noise. 

J1632-4621 

2011-06-26-07:39:22/09 

0.21 

1.709154 

562.9 

29.8 

15.2 

Relatively long period pulsar, confusion with red noise. 


J1644-44 

2011-07-03-10:59:47/11 

0.047 

0.173911 

535.1 

30.1 13.1 

Published declination not well-determined. 

Knisoel et al.|2013 

) 

J1644-46 

2011-05-17-13:02:58/02 

0.044 

0.250941 

405.8 

48.3 8.3 

Published declination not well-determined. 

Knisnel et alJ2013 

) 


J1749-2347 

2013-02-02-01:44:13/08 

0.067 

0.874486 

344.0 

28.3 

9.3 

J1809-2004 

2013-04-02-16:32:16/05 

0.18 

0.434811 

867.1 

17.9 

9.3 

J1818-1556 

2013-04-07-16:12:18/03 

0.15 

0.952709 

230.0 

30.3 

8.5 

J1746-2850 

2011-06-26-10:07:53/08 

0.12 

1.077101 

962.7 

32.2 

- 

J1644-4657 

2011-05-17-13:02:58/08 

0.098 

0.125962 

718.0 

14.7 

— 


Table C3. Previously-known binary pulsars re-detected in the 50 per cent processed data of the HTRU Galactic plane survey, sorted by their respective catalogue 
orbital periods (Porb,cat) descending order. We list the file ID (which stands for the date and the time in UT) of the relevant observation where these binaries 
were detected, as well as the detected spin period (Pobs) and the DM of the pulsar. We compare the detected S/N from the ‘standard’ (i.e., no acceleration search) 
pipeline with those from the full-length, half-length, quartered-length and one-eighth length segments of the ‘partially-coherent segmented acceleration search’ 
pipeline described in Section [3.2.31 The corresponding detected orbital acceleration (aorb) from these acceleration searches are also listed. Note that in some 
cases, only one of the two pipelines has been employed, hence the non-applicable columns are denoted by ‘=’. If the data were processed but the pulsar was not 
detected, it is denoted by ‘ x ’. 


PSR name 

-^orbjCat 

(hr) 

File ID 
(Pointing/beam) 

-^obs 

(ms) 

I^^obs 

(cm~^pc) 

‘Standard’ 

S/N 

Full-length 

S/N Uorb 

(ms-2) 

Half-length 

S/N ^orb 

(ms-2) 

V4-length 

S/N Oorb 

(ms-2) 

l/8-length 

S/N ^orb 

(ms”^) 

B1259-63 

29681.4 

2011-12-05-18:51:30/04 

47.763 

147.2 

63.5 

62.7 

0.0 

58.9 

0.03 

15.3 

202.1 

48.0 

—6.6 

J1711-4322 

22139.3 

2011-04-24-13:15:17/08 

102.618 

191.9 

18.0 

17.6 

0.06 

11.8 

-1.0 

X 

X 

X 

X 

B1800-27 

9762.7 

2011-08-16-12:50:12/12 

334.412 

162.1 

==== 

48.1 

0.0 

42.3 

1.0 

24.1 

206 

28.7 

104 

B1820-11 

8586.3 

2012-08-07-08:52:51/12 

279.837 

425.6 

==== 

47.9 

0.9 

44.7 

-2.1 

26.2 

202.1 

34.0 

44.5 

J1708-3506 

3579.2 

2011-07-07-09:57:10/03 

4.505 

146.7 

29.6 









J1750-2536 

411.4 

2013-02-06-20:15:56/10 

34.745 

179.2 

==== 

11.5 

0.064 

X 

X 

X 

X 

X 

X 

J1810-2005 

360.3 

2013-04-08-20:07:00/12 

32.822 

238.8 

9.1 

28.3 

0.33 

25.4 

-1.0 

X 

X 

18.2 

—6.6 

J1753-2240 

327.3 

2011-05-05-18:50:26/02 

95.135 

155.3 

18.5 









J1454-5846 

298.2 

2013-01-02-19:24:06/05 

45.249 

116.8 

==== 

25.5 

0.06 

19.9 

1.0 

X 

X 

X 

X 

J1543-5149 

193.5 

2011-07-04-07:23:31/07 

2.057 

51.1 

14.8 

12.9 

-0.2 

9.7 

0.03 

X 

X 

X 

X 

J1232-6501 

44.7 

2011-04-23-08:08:01/01 

88.282 

240.0 

36.5 

39.3 

0.9 

31.6 

1.0 

8.9 

308.5 

17.4 

27.4 

J1525-5545 

23.8 

2011-04-19-16:41:28/11 

11.360 

126.8 

X 

24.0 

-1.0 

32.4 

-1.0 

X 

X 

17.7 

—6.6 

B1913-l-16t 

7.8 

2010-12-30-03:53:30/01 

59.013 

169.1 

38.7 

39.7 

-1.0 

36.5 

-2.1 

9.2 

-224.2 

21.7 

—6.6 

J0737-3039At 2.5 

2010-12-30-17:11:35/01 

22.674 

47.9 

14.1 

14.1 

0.0 

32.0 

-198.9 

38.2 

-207.1 

38.0 

163.6 

B1744-24A 

1.8 

2013-02-02-01:44:13/06 

11.562 

241.9 

==== 

X 

X 

14.3 

31.9 

X 

X 

X 

X 

J1807-2459A 1.7 

2011-05-08-15:24:15/04 

3.059 

134.3 

X 

X 

X 

19.2 

-3.2 

X 

X 

14.2 

—6.6 


^ PSRs B1913+16 and J0737—3039A are technically not within the HTRU Galactic plane survey region. However, these two binaries ai'e highly relativistic and 
were observed as test pulsars at the beginning of the survey. The detected pai'ameters are listed in this table, but they are not counted towai'ds the 435 known 
pulsars re-detections. 






























